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SECTION  1 


INTRODUCTION 


This  volume  of  the  final  technical  report  contains  the 
descriptions  of  the  computer  programs  and  subprograms  which 
constitute  the  Radar  System  Simulation  Model. 

The  simulation  model  computer  program  is  divided  into 
two  distinct  phases:  (1)  Data  initialization  acivity,  and 
(2)  Simulation  activity.  In  terms  of  software,  each  phase 
is  composed  of  one  main  and  a  group  of  subprograms. 

The  subprograms  are  divided  into  five  distinct  groups: 
Stimulus /transfer  function  modules,  connection  modules, 
peripheral  modules,  supervisory  modules  and  subordinate 
modules.  These  groups  are  more  fully  described  in  the 
sections  devoted  to  each  group. 

The  flow  charts  and  cross  reference  tables  for  the 
entire  Radar  System  Simulation  Model  are  located  in  Part  2 
of  this  volume . 

The  module  descriptions,  program  listings,  flow  charts 
and  cross  reference  tables  have  all  been  cross  indexed  and 
are  located  in  Section  10.  This  section  has  been  included 
in  both  parts  of  the  volume  for  convenience. 
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SECTION  2 


SIMULATION  DATA  LOADER  EXECUTIVE 

(M  A  I  N  -  1) 

The  simulation  data  loader  serves  as  the  interface  be¬ 
tween  the  user  and  the  segment  of  the  computer  program  which 
performs  the  simulation.  This  segmentation  was  necessary 
since  input  data  in  a  format  convenient  to  the  user  must  be 
converted  to  a  form  suitable  for  use  by  the  simulation  mod¬ 
ules.  The  punch  cards  which  define  the  simulation  to  be 
performed  are  of  two  types:  simulation  control  cards  and 
module  parameter  data  cards.  The  simulation  control  cards 
determine  what  operations  are  to  be  performed  in  the  simula¬ 
tion  activity.  This  includes  not  only  the  scheduling  of  mod¬ 
ules  for  execution  but  also  the  movement  of  data  to  and  from 
temporary  storage,  and  the  modification  of  parameters  for 
multiple  executions  of  a  simulation  model  configuration. 

Each  control  card  with  the  exception  of  ENDPAS ,  ENDCFG,  and 
ENDJ9B  is  converted  into  a  control  word  by  the  control  word 
generator  and  placed  in  the  control  word  block.  Each  control 
word  contains  the  code  number  of  the  operation  to  be  performed, 
the  module  to  be  executed  and  data  set  reference  number  if 
required.  The  module  parameter  data  cards  define  the  para¬ 
meters  of  the  simulation  modules  to  be  executed  in  the  simu¬ 
lation.  Each  module  requiring  input  data  has  a  unique  name- 
list  which  contains  the  input  parameters. 

Figure  2-1  is  a  block  diagram  of  the  Initializer  load 
module.  In  the  block  diagram  the  data  flow  paths  are  shown 
as  solid  lines  and  control  paths  are  shown  as  dashed  lines. 
Control  normally  is  retained  by  the  Simulation  Data  Loader 
Executive  but  control  transfer  does  occur  to  the  subprograms 
CLINT  and  PHENC .  The  blocks  containing  the  letter  M  represent 
data  transfers  between  two  storage  areas.  Arrows  are  used  to 
indicate  the  direction  of  data  flow.  The  blocks  containing 
C&M  perform  perform  a  search  operation  to  determine  if  the 
word  stored  in  a  buffer  is  a  member  of  a  reference  dictionary. 
If  the  search  is  successful,  the  data  in  the  buffer  is  trans¬ 
ferred  to  the  storage  area  designated  by  an  arrow.  The  func¬ 
tion  served  by  certain  blocks  in  the  diagram  are  evident  from 
their  titles.  Those  blocks  requiring  further  explanation 
are  discussed  in  the  following  paragraphs. 
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RECORD  NUMBER 


SIMULATION  DATA  LOADER  EXECUTIVE 


The  ModuLe  Parameter  Namelist  Dictionary  serves  the 
function  of  directing  the  input  data  to  the  proper  location 
within  the  Module /System  Parameter  Table.  Appendix  B  con¬ 
tains  a  list  of  the  parameters  in  each  namelist  and  the 
storage  location  assigned  to  each  parameter.  The  namelist 
name  is  the  same  as  the  Module  Reference  Number,  i.e.,  the 
namelist  for  module  101  is  NL101. 

The  Module  Reference  Number  Dictionary  is  a  cross 
reference  between  Module  Reference  Numbers  and  the  overhead 
operations  which  must  be  accomplished  prior  to  execution  of 
a  module.  For  example,  to  generate  a  phase  encoded  waveform 
requires  that  the  user  supplied  data  be  preprocessed  by 
sub  rou  ti ne  PHENC . 

The  Command  Word  Dictionary  is  a  cross  reference  between 
input  command  words  and  the  operations  to  be  performed  in 
the  simulation.  The  control  word  is  entered  in  columns  1 
through  6  of  the  simulation  control  card.  The  following  is 
a  list  of  the  command  words  and  the  corresponding  operation 
initiated  by  each: 

EXEC  This  command  word  schedules  execution  of  the 

module  corresponding  to  the  number  appearing  in 
columns  15  through  17  of  the  control  card.  No 
input  data  cards  are  required. 

LDEXEC  This  command  word  is  the  same  as  EXEC  with  the 
additional  requirement  that  input  data  is  read. 

The  namelist  name  for  entering  the  data  is  the 
same  as  the  module  reference  number  contained  in 
columns  15  through  17  of  the  simulation  control 
card . 

M9DIFY  This  command  word  causes  the  data  loaded  in  a 

previous  step  to  be  modified.  This  is  typically 
used  to  change  parameters  when  multiple  simula¬ 
tion  passes  are  to  be  made.  The  step  number  to 
be  changed  is  entered  in  columns  10  through  12  of 
the  simulation  control  card,  right  adjusted.  The 
Module  Reference  Number  is  entered  in  columns 
15  through  17  of  the  simulation  control  card. 

ENDPAS  This  command  word  signifies  the  end  of  a  pass 
through  a  simulation  configuration. 
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ENDCFG  This  command  word  signifies  the  end  of  a 
configuration. 

ENDJ9B  This  command  word  signifies  the  end  of  a 
simulation  job. 

S T9 REX  This  command  word  causes  the  contents  of  the 

XT  signal  storage  array  to  be  stored  on  a  tempo¬ 
rary  data  set,  usually  a  disc  file.  The  Data 
Set  Reference  Number  is  entered  in  columns  11 
and  12  of  the  simulation  control  card. 

ST9REY  Same  as  ST9REX  except  the  contents  of  signal 
storage  array  YT  are  stored. 

ST9REA  Same  as  ST9REX  except  the  contents  of  auxiliary 
storage  array  XA  are  stored. 

ST9REB  Same  as  STOREX  except  the  contents  of  auxiliary 
storage  array  XB  are  stored. 

L9AD  X  This  command  word  causes  the  data  located  in  a 

temporary  file  to  be  loaded  into  the  signal  stor¬ 
age  array  XT.  The  Data  Set  Reference  Number  is 
entered  in  columns  11  and  12  of  the  simulation 
control  card. 

L9AD  Y  Same  as  L9AD  X  except  the  signal  storage  array 
YT  is  loaded. 

L9AD  A  Same  as  L9AD  X  except  the  auxiliary  storage 
array  XA  is  loaded. 

L9AD  B  Same  as  L9AD  X  except  the  auxiliary  storage 
array  XB  is  loaded. 

REWIND  This  command  word  causes  a  temporary  file  to  be 
rewound.  For  a  disc  file  this  moves  the  access 
pointer  to  the  beginning  of  the  data  set.  The 
Data  Set  Reference  Number  is  entered  in  columns 
11  and  12  of  the  simulation  control  card. 

The  output  from  the  data  loader  are  stored  on  temporary 
storage  files.  Data  Set  Reference  Number  11  contains  the 
control  word  block  which  controls  the  simulation  activity 
and  the  module/system  parameter  tables  used  to  define  the 
modules  and  system  characteristics.  Data  Set  Reference  Number 
2  contains  the  clutter  scatterer  parameters. 
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Though  reduction  is  defined  as  a  post  processing  operation, 
it  is,  in  reality,  more  of  a  resolution  operation  and, 
therefore,  should  be  executed  within  the  resolution  pro¬ 
cedure  or  prior  to  other  post  processing  operations. 


(6)  Ordering 

Literals  to  the  right  of  the  right-most  framed  literal  are 
free  to  be  reordered.  Reordering  may  not  occur  across 
framed  literals  without  special  care. 


2.4  FACT  DETERMINATION 


OL-deduction ,  as  has  been  stated,  is  a  refutation  procedure.  A 
questionable  assertion,  that  is,  a  query,  is  presented  in  negated  form  to 
the  OL-deduction  mechanism.  Deductions  are  then  generated  and  are  sought 
to  be  refuted.  A  refutable  deduction  implies  the  existance  of  a  collection 
of  rules  and  facts  (possibly  a  null  collection)  that  refutes  the  negated  query, 
and  consequently,  satisfies  the  assertion  of  the  query.  It  is  the  collection 
of  facts  that  is  the  answer  to  the  positive  query  and  that  is  provided  by  the 
inference  system, 

A  resolved  deduction  is  a  clause  consisting  of  only  resolution 
literals;  that  is,  it  is  a  deduction  clause  that  has  been  fully  resolved. 

Each  resolution  (framed)  literal  in  a  resolved  deduction  clause 
is  a  literal  about  which  the  information  system  can  supply  facts  or  is  an 
acceptable  inferred  literal.  An  inferred  literal  will  often  be  a  simple 
binary  proposition,  signifying  a  yes  or  no;  true  or  false;  0  or  1  value.  In 
the  more  general  case,  where  the  framed  literal  is  a  complex  predicate,  the 
literal  will  tend  to  be  a  fact  literal  for  which  the  information  system  will 
supply  all  facts.  For  example,  given  |  L(X,5)  I  where  L(X,5)  means  "X"  is  a 
worker  in  factory  "5",  then  the  information  system  will  search  out  all  workers 
in  factory  number  5  and  supply  these  as  facts  satisfying  the  literal.  We  call 
this  procedure  "fact  determination”  When  fact  determination  is  applied  to 
all  of  the  fact  literals  of  a  completely  resolved  deduction,  we  then  determine 
a  largest  subset  of  these  facts  that  can  satisfy  the  deduction  as  an  entirety. 
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SECTION  3 


SIMULATION  CONTROLLER 
(M  A  I  N  -  2) 

The  simulation  controller  serves  as  the  switchboard 
which  connects  the  simulation  modules  together  in  the  manner 
prescribed  by  the  control  cards  read  by  the  simulation  data 
loader.  Figure  3-1  is  a  block  diagram  of  the  Simulation 
Load  Module.  The  control  word  blocks  and  System/Module  Para¬ 
meter  Tables  are  stored  in  DSRN  #11  which  was  initialized 
during  the  simulation  data  loader  activity.  The  clutter 
scatterer  parameters  are  stored  in  DSRN  #2  which  was  initial¬ 
ized  during  the  simulation  data  loader  activity.  In  addition, 
provisions  are  available  for  reading  and  writing  data  on 
user  defined  disc  and  magnetic  tape  data  sets. 

There  are  two  storage  arrays  allocated  for  storing  sim¬ 
ulated  signal  data.  These  arrays  are  designated  as  XT  and 
YT.  In  addition,  two  auxiliary  arrays  are  allocated  for 
storing  either  signal  data,  processed  data  such  as  probability 
distributions,  or  other  miscellaneous  data.  These  arrays  are 
designated  as  XA  and  XB.  All  four  arrays  are  composed  of 
8192  elements. 

The  following  is  a  description  of  the  sequence  of  events 
that  occur  in  performing  a  step  of  the  simulation. 

1.  The  operation  to  be  performed  in  the  Jth  step  is 
defined  by  the  control  word  stored  in  location  J 
of  the  control  word  block.  This  word  is  com¬ 
pared  with  the  control  word  dictionary  to  deter¬ 
mine  (1)  the  operation  to  be  performed,  (2)  the 
reference  number  of  the  module  to  be  executed, 
and  (3)  the  Data  Set  Reference  Number. 

2.  If  the  operation  to  be  performed  is  the  execution 
of  a  module  with  no  input  data  requirements,  then 
control  is  transferred  to  the  designated  module. 

If  the  operation  to  be  performed  is  execution  of 
a  module  requiring  input  data,  then  the  System/ 
Module  Parameter  table  is  loaded  from  DSRN  #11 
and  control  is  transferred  to  the  designated 
module.  For  operations  involving  data  transfers 
the  contents  of  the  designated  array  are  either 
loaded  from  or  transferred  to  the  DSRN  which  was 
specified . 
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Figure  3-1  SIMULATION  LOAD  MODULE  FUNCTIONAL  BLOCK  DIAGRAM 


3.  Once  the  operation  for  the  Jth  step  is  completed 
the  job  step  counter  is  incremented  by  1  and 
the  process  is  repeated. 

When  a  blank  (all  zeros)  control  word  is  encountered, 
the  simulation  pass  is  ended.  The  job  step  counter  is  reset 
to  1  and  a  new  pass  through  the  simulation  is  begun  with  a 
different  set  of  parameters.  A  block  of  data  containing  all 
zeros  is  used  to  designate  the  end  of  a  simulation  configura¬ 
tion.  When  this  is  encountered  the  simulation  controller 
will  move  the  next  block  of  data  into  the  control  word  buffer. 
This  new  block  of  control  words  defines  a  new  simulation  con¬ 
figuration.  When  a  new  control  word  block  containing  all 
zeros  is  encountered,  the  simulation  job  is  terminated. 

Flow  Chart:  Page  9-34 

Cross  Reference  Table:  Page  9-210 
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CALL  DTOA  (IX,X) 

Where:  IX  contains  the  Input  Waveform 

X  contains  the  Output  Waveform 

5.  RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  If  the  computed  digital  output  is  greater  than  the 
largest  number  which  can  be  specified  by  NBITS,  it 
will  be  set  to  the  maximum  number. 

b.  Flow  Chart:  Page  9-98 

c.  Cross  Reference  Table:  9-220. 

6.  THEORY  OF  OPERATION 


The  analog- to-digital  conversion  is  performed  by  divid¬ 
ing  the  floating  point  input  signal  by  the  value  of  the  least 
significant  digital  bit.  If  round  off  is  to  be  performed, 

0.5  is  then  added  to  the  value.  The  value  thus  obtained  is 
converted  to  an  integer  number  and  its  absolute  value  is 
tested  against  the  maximum  allowable  number,  2**(NBITS-1) . 

If  the  magnitude  of  the  signal  is  greater  than  the  maximum 
allowable,  it  is  set  equal  to  the  maximu.  The  basic  mechani¬ 
zation  equation  for  the  ATOD  module  is: 

IX  (J)  =  IFIX(X( J) /XLSB  +  0.5  *  IROFF) 

The  digital- to-analog  conversion  is  made  by  changing 
the  input  to  a  floating  point  number  and  multiplying  it 
times  the  value  of  the  least  significant  bit.  The  basic 
mechanization  equation  for  the  DTOA  module  is: 

X(J)  =  FLOAT  (IX(J))  *  ZLSB 
where:  ZLSB  =  BOOL  (IX(N196)) 
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SUBROUTINE  CDIGFL 


1.  MODULE  IDENTIFICATION: 


Name 

Classification  Code 

Reference  Number 

CDIGFL 

LTI-4 

403 

CDFNCL 

LTI-4 

404 

2.  PURPOSE: 


This  subroutine  is  used  to  simulate  a  cross-coupled 
digital  filter  and  synthesize  a  filter  transfer  function 
having  one  complex  pole  and/or  one  complex  zero. 


3.  INPUT  PARAMETERS: 
Name  0/R  T 

FFR  0  F 

FFI  0  F 

FBR  0  F 

FBI  0  F 


Description 

Feed  -  Forward  coefficient  -  R 
Feed  -  Forward  coefficient  -  I 
Feedback  coefficient  -  R 
Feedback  coefficient  -  I 


4.  CALLING  SEQUENCES: 


CALL  CDIGFL  (X,Y) 


Where:  X  contains 

Y  contains 
X  contains 

Y  contains 


the  Input  Waveform  -  R 
the  Input  Waveform  -  I 
the  Output  Waveform  -  R 
the  Output  Waveform  -  I 


The  storage  registers  (TX1  and  TY1)  are 
cleared  before  execution  begins. 


CALL  CDFNCL  (X,Y) 


Where :  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 
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The  storage  registers,  TX1  and  TY1,  are  not 
cleared  before  execution  begins. 


I 

I 


l 


5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-96 

b.  Cross  Reference  Table:  9-220 

6.  THEORY  OF  OPERATION 


The  block  diagram  of  the  complex  digital  filter  simu¬ 
lated  by  this  module  is  shown  in  Figure  CDIGFL-1.  The  Z- 
plane  transfer  function  is  given  by  the  following  expression: 


TCZ'i  =  7  ~  (FFR  +  iFFI) 

'  '  n  /■C'Tjn  _L  !PDT\ 


Z  -  (FBR  +  jFBI) 


FFR 
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SUBROUTINE  CFAR 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

CFAR  NL-2  459,460 

2 „  PURPOSE: 

This  subroutine  is  used  to  simulate  a  constant 
false  alarm  rate  (CFAR)  video  processor. 

3.  INPUT  PARAMETERS: 

Name  O/R  T  Description 

TAVG  R  F  Width  of  averaging  window  used 

in  determining  video  gain 
(TAVG  >  2*TI) 

4.  CALLING  SEQUENCES: 

CALL  CFAR(VIN, VOUT) 

Where:  VIN  contains  the  Input  Waveform 

VOUT  contains  the  Output  Waveform 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  input  and  output  arrays  cannot  be  equivalenced. 

b.  Flow  Chart:  Page  9-66 

c.  Cross  Reference  Table:  Page  9-216 

6 .  THEORY  OF  OPERATION 


The  block  diagram  of  the  CFAR  processor  simulated 
by  this  module  is  shown  in  Figure  CFAR-1.  Each  sample 
of  the  input  waveform  is  passed  through  a  variable 
gain  amplifier  which  is  controlled  by  the  average  value 
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BLOCK  DIAGRAM  OF  CFAR 


of  this  waveform  surrounding  the  sample  being  processed. 
The  gain  applied  to  the  Jth  sample  is  given  by  the 
following  expression: 


G(j\  -  _ 2*NCELL2 _ 

K  '  K  =  J+NCELL2 

]P(|VIN(K)|  -  |  VIN (J ) | ) 

K  =  J-NCELL2 

where:  NCELL2  <  J  <  NPTS  -  NCELL2 

NCELL2  =  IFIX(TAVG*0„ 5/TI) 

TI  =  independent  variable  increment  between 
samples 

NPTS  =  number  of  samples  in  the  input  array 
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SUBROUTINE  CLUTTR 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

CLUTTR  LTI-3  503 

or 
LTV- 3 

2.  PURPOSE : 

This  subroutine  generates  the  transfer  function  of 
a  clutter  volume.  A  pictorial  diagram  of  the  clutter  vol¬ 
ume  is  shown  in  Figure  CLUTTR-1. 

3.  INPUT  PARAMETERS: 

All  required  input  data  is  entered  v^'a  the  clutter 
model  initialization  subroutine  (CLINT). 


4. 


CALLING  SEQUENCES: 


CALL  CLUTTR (XT, YT.GAZ, GEL) 


WHERE: 


XT  contains  the  Output  Impulse 
response  -  R 


YT  contains  the  Output  Impulse 
response  -  I 


GAZ-Temporary  storage  for  antenna 
azimuth  gain  values  (MM=number  of 
storage  cells  required) 


GEL-Temporary  storage  for  antenna 
elevation  gain  values  (NN  -  number  of 
storage  cells  required) 


a.  The  subroutines  CLINT  must  have  been  successfully 
executed  prior  to  calling  CLUTTR.  See  p  6-4  for 
restrictions  on  CLINT. 

b.  The  functions  AZGAIN  and  ELGAIN  are  used  to 
compute  the  antenna  gain  associated  with 
each  scatterer. 
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MODEL  GEOMETRY 


c.  Flow  Chart:  Page  9-112 

d.  Cross  Reference  Table:  Page  9-222 

6.  THEORY  OF  OPERATION 

The  clutter  model  is  based  on  the  premise  that  a 
volume  of  clutter  can  be  represented  by  an  ensemble  of 
discrete  scatterers.  The  RCS  of  these  scatterers  is 
derived  from  a  theoretical  probability  density  function, 
but  can  be  easily  modified  to  allow  the  use  of  a  density 
function  derived  from  measurements  data.  The  phase  behavior 
of  the  scatterer  is  determined  from  motion  of  the  scatterer 
due  to  wind  and  from  any  inherent  frequency  spread  derived 
from  measurements  data.  The  basic  mechanization  equation 
for  the  clutter  model  is  given  by  the  following  expressions: 


a .  MODE  = 


XT(K)  4-  jYT(K)  = 
where:  GEL(N)  - 


I  (stationary  clutter) : 

NN  MM 

£  GEL (N)  £  GAZ(M)  CLUX(IFNT)+jCLUY(IPNT) 
N-l  M-l 

Normalized  antenna  gain  corresponding 

to  an  elevation  angle  of  EL000  +  DELEL  *(N-1) 


GAZ(M)  -  Antenna  gain  corresponding  to  an  azimuth 
angle  of  AZ000  +  DELAZ  *  (M-l) 

IFNT  -  Scatterer  element  pointer  which  is 
determined  as  follows: 

I PNT  --  (K*(N-1)  (M-1)  +  (N-1)*MM+M)  MODULO  250 

CLUX(  )-  Array  containing  real  part  of  clutter 
scatterer  specification 

CLUY(  )-  Array  containing  imaginary  part  of 
clutter  scatterer  specification 


b.  MODE  =  2 
NN 

XT (K)+jYT (K)  =  £ 
N=1 


(Time  varying  clutter): 

m  ,  ,  PHDEL 

GEL(N)  £  GAZ(M)|CLUX(PDNT)+jCLUY(IDNT)|*  e 


"New”  CLUX (IDNT)+ jCLUY (TDNT) 
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Where:  PHD EL  =  is  the  phase  term  to  be  applied  to  each 

scatter 


=  57, 


.29573*DELTIM*DOPFRQ*COSj|xVLANG-DELAZ*<M-l)]y57 .29567J  +  PHRW 


DELTIM  =  Radar  interpulse  period 

PHRW  =  The  random  walk  phase.  This  variable  is 
generated  by  a  random  number  generator 
from  a  umiform  distribution. 


-RWPH  <  PHRW  <  RWPH 

2.0  2.0 


The  pulse  phase  shift  due  to  doppler  is  computed  via 
a  recursive  loop  which  is  more  efficient  than  a  direct 
evaluation  of  the  cosine  function. 


C(N+1)  =  C(N)*DELCS  -  S(N)*DELSN 
S (N+l )  -  C(N)*DELSN+S(N)*DELCS 

where  C  =  Phase  shift  due  to  doppler 

C(l)  -  D0PFRQ*DELTIM*2 n  *  COS(XVLANG) 

S(l)  =  DOPFRQ*DELTIM*2t7  *  SIN  (XVLANG) 

The  clutter  scatterer  parameters  (CLUX  and  CLUY)  are 
updated  by  PHDEL  for  each  execution  of  the  clutter  model 
and  the  new  values  are  stored  in  the  random  disc  file. 
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SUBROUTINE  DCFAR 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

DCFAR  NL-2  440,441 

2.  PURPOSE: 

This  subroutine  is  used  to  simulate  a  digital  constant 
false  alarm  rate  (CFAR)  video  processor. 

3.  INPUT  PARAMETERS: 

Name  O/R  T  Description 

NCELL  R  I  Number  of  cells  to  be  averaged 

in  computing  the  video  gain. 

4.  CALLING  SEQUENCES: 

CALL  DCFAR  (IN,  IOUT) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  input  and  output  arrays  cannot  be  equivalenced. 

b.  All  calculations  are  in  integer  format. 

c.  Flow  Chart:  Page  9-93 

d.  Cross  Reference  Table:  Page  9-220 

6 .  THEORY  OF  OPERATION 


The  block  diagram  of  the  digital  CFAR  processor  simu¬ 
lated  by  this  module  is  shown  in  Figure  DCFAR- 1.  Each 
sample  of  the  input  waveform  is  passed  through  a  variable 
gain  amplifier  which  is  contro  led  by  the  average  value  of 
the  waveform  surrounding  the  sample  being  processed.  The 
gain  applied  to  the  Jth  sample  is  given  by  the  following 
expression: 
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_ NCELL _ 

IG  =  K=J+NC£LL2 

£  I  IN (K)|  -  IN(J)| 

K=J-NCELL2 

where:  NCELL  2  <  J  <  NPTS  -  NCELL2 

NCELL  2  =  NCELL/ 2 

NPTS  =  Number  of  input  samples  in  the  input 
array 

The  value  of  the  Jth  output  sample  is  computed  as  follows 

IOUT(J)  =  IN ( J )  *  IG 
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SUBROUTINE  DFT 


1.  MODULE  IDENTIFICATION: 


Name 


Classification  Code  Reference  Number 


DFT 

DFTF0 

DFTRF 

DFTNCL 


LTI-2 

201 

LTI-2 

234 

LTI-2 

233 

LTI-2S 

None 

2.  PURPOSE: 


This  subroutine  computes  the  finite  discrete  Fourier 
transform  of  a  sequence  of  8 -functions. 


3.  INPUT  PARAMETERS: 


Name  O/R  T 

SIMF0  0  F 

FI  R  F 

FEXT  R  F 


Description 

Center  frequency  of  discrete  output 
spectrum  (DFTF0  only) 

Frequency  increment  between  samples 
of  the  output  spectrum 

Width  of  output  spectrum 


NIMP  R  I 


Number  of  8-functions  to  be  trans¬ 
formed 


IDMP  R  I 

F0  OF 

I NORM  R  I 


Diagnostic  data  dump  parameter 
=  0;  no  diagnostic  data 
■  1;  dump  internal  parameters:  PS,  DELPS, 
NPTS,  NIMP,  DIN 

«  2;  dump  output  spectrum:  PS,  DELPS, 
NPTS,  NIMP,  DIN 

Center  frequency  of  discrete  output 
spectrum  (DFTRF  only) 

Normalization  Flag 

-  0;  fN  -  TI 

=  1;  fN  -  1.0/NIMp 

-  2;  fN  -  TI 

*  3;  fN  »  1’0/NIMP 

-  4;  fN  -  1.0 
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Name  0/R  T 


Description 


TI  OF  Time  increment  between  simulated  wave¬ 

form  samples.  Used  in  this  subroutine 
for  normalization  only  (required  only 
if  INORM  =  0  or  2) 


DIN  R  F  Array  containing  the  parameters  of 

the  8 -functions.  The  specification 
for  the  Jth  8-function  is  the  following: 
DIN(1,J)  =  phase  angle 
DIN(2,J)  =  time  of  occurrence 
DIN(3,J)  =  amplitude 


4.  CALLING  SEQUENCES: 

Video  spectrum 
CALL  DFT  (X,Y) 

Where:  X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  range  of  the  independent  variable,  f,  used  in 
computing  the  spectrum  is: 

FEXT  <f  FEXT 
"2  2 


IF  spectrum 
CALL  DFTF0  (X,Y) 

Where:  X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  range  of  the  independent  variable,  f,  used  in 
computing  the  spectrum  is: 

SIMF0  -  <  f  <  SIMF0  + 


RF  spectrum 
CALL  DFTRF  (X,Y) 
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where:  X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  range  of  the  independent  variable,  f,  used  in 
computing  the  spectrum  is: 

F0  -  SSI  <  f  <  F0  +  E|XT 

Transform  additional  sample 
CALL  DFTNCL  (X,Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  computed  spectrum  is  added  to  the  spectrum  contained 
in  arrays  X  and  Y.  The  range  of  the  independent  variable 
is  determined  by  the  previous  execution  of  DFT,  DFTF0, 
or  DFTRF. 


RESTRICTIONS. 


JUIREMENTS,  MISCELLANEOUS  DATA 


a.  This  subroutine  is  quite  time  consuming  and  therefore 
should  be  used  only  when  necessary.  For  most  appli¬ 
cations  the  AFFT  subroutine  which  uses  the  Fast 
Fourier  Transform  algorithm  should  be  used.  The 
execution  time  for  this  subroutine  can  be  estimated 
using  the  following  expression 


FEXT 

FI 


x 


NIMP  x  50 


b.  This  subroutine  is  used  normally  in  conjunction  with 
other  subroutines  such  as  PHDEC  and  CGEN. 


c.  The  output  spectrum  has  a  center  frequency  of  zero, 
i.e.  it  is  the  low  pass  equivalent. 

d.  Flow  Chart:  Page  9-142 

e.  Cross  Reference  Table:  Page  9-226 

f.  DFT  lets  user  generate  100  delays  to  be  considered. 
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6. 


theory  OF  operation 


expression: 


* 

f-  FSB?  a,  s  - j2uft 

S(£)-F  J  C 

m=l 

^  ... _  not 


where :  FN  is  the  normalization  factor  determined  by  the 

“tfrcSng  the  order  of  integration  and  auction 
the  following  result  xs  obtained. 


N1MP 

s(£)  =  %  Z-r 

m=l 


j0(m)  - j2»ft(m) 

a/^,'»e  e 


If  only  K  equally  spaced  samples  of  S(f)  are  computed, 
the  following  expression  results. 

-  fY  a(m)  e3*«  e-j2.mt(m) 

Sk  -  s*(f)  -  F»Z-,  aW 


>f=kFI+f0 


tn  =*1 


K  _  FEXT  d  f  ^  defined  in  the  following 
where:  K  "  FI  7 ,  ° 
manner  for  each  module  name. 


FEXT 

DFT:  fD  -  ‘  ~T~ 


.  fext 

DFTF0:  fQ  ~  SIMF0  ~  2 


DFTRF :  f0  =  F0  - 2 


FEXT 
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SUBROUTINE  DIGFIL 


1.  MODULE  IDENTIFICATION: 


Name 

Classification  Code 

Reference  Number 

DIGFIL 

LTI-4 

461 

DIGFNC 

LTI-4 

462 

ECMFL 

LTI-4 

None 

2.  PURPOSE: 

This  subroutine  is  used  to  simulate  a  multiple  section 
digital  filter. 


3 .  INPUT  PARAMETERS : 
Name  O/R  T 

NP  R  I 

SF  R  F 

FB  OF 


FF  OF 


Description 

Number  of  2-delay  sections 
Scale  factor 

Array  containing  the  feedback 
coefficients.  The  coefficients 
for  the  Kth  section  are  the 
following: 

FB(1,K)  =  1-delay  feedback  coef. 
FB(2,K)  =  2-delay  feedback  coef. 

Array  containing  the  feed-forward 
coefficients.  Tine  coefficients 
for  the  Kth  section  are  the 
following: 

FF(1 ,K)  =  0-delay  feed-forward 
coef. 

FF( 2 , K)  =  1-delay  feed- forward 
coe  f . 


4 .  CALLING  SEQUENCES ; 
CALL  DIGFIL  (X,Y) 
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Where : 


X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  storage  register  arrays  (XM  and  YM) 
are  cleared  before  execution  begins. 

CALL  DIGFNC  (X,Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  storage  register  arrays  (XM  and  YM) 
are  not  cleared  before  execution  begins. 

CALL  ECMFL  (X,Y) 

Where:  X  contains  the  Input  Waveform  sample  -  R 

Y  contains  the  Input  Waveform  sample  -  I 

X  contains  the  Output  Waveform  sample  -  R 

Y  contains  the  Output  Waveform  sample  -  I 

One  complex  sample  processed  for  each 
execution  of  the  module . 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-101 

b.  Cross  Reference  Table:  Page  9-221 

6.  THEORY  OF  OPERATION 


The  block  diagram  of  one  section  of  the  digital  filter 
simulated  by  this  module  is  shown  in  Figure  DIGFIL-1.  The 
Z-plane  transfer  function  for  this  section  is  given  by  the 
following  expression: 


Tk(Z)  =  FF(1,K) 


22  +  z  + 

Z2  -  FB(1 ,K)  Z  -  FB(2,K) 


The  Z-plane  transfer  function  for  the  complete  filter  is 
given  by  the  following  expression: 

K>NP 

T(Z)  =  SF  (I  Tk(Z) 

K-l 


FF(1,K) 


Figure  DIGFIL-1  BLOCK  DIAGRAM  OF  DIGFIL/DIGFNC 
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SUBROUTINE  DIGFSF 


1 .  MODULE  IDENTIFI CATION : 


Name 


Classification  Code  Reference  Number 


DIGFSF  LTI-4  463 

ECMFSU  LTI-4  None 

ECMFSF  LTI-4  None 

2.  PURPOSE : 

This  subroutine  is  used  to  simulate  a  digital  filter 
based  on  the  frequency  sampling  design  concept. 


3.  INPUT  PARAMETERS: 
Name  0/R  T 

RADIUS  R  F 

NSAM  R  I 

NDFZ  R  I 

FSAM  R  F 


Description 

Radius  of  circle  in  Z-plane  which 
is  locus  of  pole/zero  locations 

Number  of  frequency  samples  used 
to  specify  the  transfer  function 

Number  of  zeros  synthesized 

Array  containing  the  filter 
transfer  function  specification. 
The  coefficients  for  the  Kth 
sample  are  the  following: 

FSAM  (1,K)  =  real  coefficient 
FSAM  (2,K)  =  imaginary  coefficient 

The  frequency  associated  with  the 
Kth  sample  is  given  by  the  follow¬ 
ing  expression. 

fK  *  nSfz  <k-NDC>  K  "  1.  NSAM 

where:  NDC  -  IFIX  (NSAM/2)  +  1 

fg  *  Sampling  rate  of 

data  to  be  processed 
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CALLING  SEQUENCES: 


Process  a  waveform 
CALL  DIGFSF  (X,Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

The  storage  register  arrays  (S  and  PF) 
are  cleared  before  execution  begins. 


Set  up  for  use  by  ECM  module 
CALL  ECMFSU  (X,Y) 

Where:  X  contains  the  Input  Waveform  sample  - 

Y  contains  the  Input  Waveform  sample  - 
X  contains  the  Output  Waveform  sample  - 

Y  contains  the  Output  Waveform  sample  - 

The  storage  register  arrays  (S  and  PF) 
are  cleared  and  coefficients  are 
ini tialized . 

Process  1  sample  (called  by  ECM  module) 

CALL  ECMFSF  (X,Y) 

Where:  X  contains  the  Input  Waveform  sample  - 

Y  contains  the  Input  Waveform  sample  - 
X  contains  the  Output  Waveform  sample  - 

Y  contains  the  Output  Waveform  sample  - 

The  storage  register  arrays  (S  and  PF) 
are  not  cleared  before  execution  begins 

RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 


a.  This  subroutine  is  structured  to  simulate  a  low 
pass  digital  filter  only. 
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b.  References: 


C.  M.  Rader  and  B.  Gold:  "Digital  Filter  Design 
Techniques  in  the  Frequency  Domain",  Proc.  IEEE. 
Vol .  55,  pp  149-171,  Feb.  1967. 

L.  R.  Rabiner,  B.  Gold,  and  C.  A.  McGonegal: 

"An  Approach  to  the  Approximation  Problem  for 
Nonrecursive  Digital  Filters",  IEEE  Trans.  Audio 
Electroacoust. .  Vol.  AU-18,  pp  83-106,  June  1970. 

L.  R.  Rabiner  and  B.  Gold:  Theory  and  Application 
of  Digital  Signal  Processing,  Prentice-Hall,  Inc., 
Englewood  Cliffs,  N.  J.,  1975,  pp  105-123. 

c.  Flow  Chart:  Page  9-107 

d.  Cross  Reference  Table:  Page  9-222. 

6.  THEORY  OF  OPERATION 


The  basic  concept  of  this  filter  design  technique  is 
as  follows.  A  comb  filter  is  used  to  synthesize  N  zeros 
around  a  circle  of  radius  R  in  the  Z-plane.  The  banks  of 
parallel  resonators  are  placed  in  cascade  with  the  comb 
filter.  The  resonators  are  chosen  such  that  once  the  des¬ 
ired  passband  region  their  poles  exactly  cancel  the  zeros. 
The  weight  applied  to  the  resonators  are  samples  of  the 
desired  frequency  domain  transfer  function. 

The  block  diagram  of  the  digital  filter  simulated  by 
the  subroutine  is  shown  in  Figure  DIGFSF-1.  Since  the 
phases  at  resonance  of  the  consecutive  lesonators  differ 
by  ,  a  phase  adjustment  is  inserted  after  each  resonator. 
The  location  of  the  zeros  synthesized  by  the  comb  filter 
are  shown  in  Figure  D1GFSF-2.  The  block  diagram  of  a  typi¬ 
cal  complex  resonator  is  shown  in  Figure  DIGFSF-3. 
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Figure  DIGFSF-1  BLOCK  DIAGRAM  OF  DIGFSF 
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SUBROUTINE  DIGTFL 


1„  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

DIGTFL  NL-2  422,423 

2  -  PURPOSE: 

This  subroutine  is  used  to  simulate  a  digital 
transversal  filter,, 

3.  INPUT  PARAMETERS : 


Name 

O/R 

T 

Description 

NTAPS 

R 

I 

Number  of  taps 

ITAP 

R 

I 

Array  containing  the  tap  specifi¬ 
cations  „  The  coefficient  for  the 
Kth  tap  are  the  following: 

ITAP(1,K)  =  numerator  of  tap  gain 

ITAP(2,K)  =  denominator  of  tap  gain 

ITAP(3,K)  =  delay  in  sampling 
increments 

4 •  CALLING  SEQUENCES : 

CALL  DIGTFL (IN, IOUT) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 

5 ,  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a-  This  device  is  typically  used  as  a  digital  phase 
decoder 

b„  Flow  Chart:  Page  9-95 

Co  Cross  Reference  Table:  Page  9-220 
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6 .  THEORY  OF  OPERAT ION 


The  block  diagram  of  the  digital  transversal  filter 
simulated  by  this  module  is  shown  in  Figure  DIGTFL-1. 

All  calculations  performed  by  this  module  are  in  integer 
format. 
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DIGTFL-1  BLOCK  DIAGRAM  OF  DIGTFI. 
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SUBROUTINE  ECM 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

ECM  S0U-1S  512 


2.  PURPOSE: 

This  subroutine  simulates  a  noise  modulation  jammer. 


3 .  INPUT  PARAMETERS : 


Name 

TIME 

JRNG0 

JRSIM 

JMAZ 

JHGT 

JERP 

JFMBW 

JPW 

JF0 

JVEL 

JPER0D 


0/R  T 
R  F 

0  F 

R  F 

R  F 
0  F 
R  F 
0  F 
0  F 

0  F 

0  F 

0  F 


Description 

Elapsed  time  since  beginning  of 
simulation 

True  jammer  range  corresponding 
to  TIME  =  0.0 

Jammer  range  to  be  used  in  com¬ 
puting  the  starting  time  of  the 
j  ammer . 

Jammer  azimuth  angle 

Jammer  height  above  ground 

Jammer  effective  radiated  power 

Jammer  swept  bandwidth 

Jammer  on  time  for  pulsed 
j  ammer s 

Jammer  center  frequency  offset 
from  radar  center  frequency 

Jammer  radial  velocity  with 
respect  to  the  radar  site 

Jammer  pulse  repetition  interval 
for  pulsed  jammer 
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Name 


N2 

NDFZ 


0/R  T  Description 

R  I  2**N2  =  total  number  of  signal 

storage  array  elements  available 

R  I  Number  of  zeros  to  be  synthe¬ 

sized  by  the  frequency  sampling 
digital  filter  (DIGFSF) . 


In  addition  to  the  above  parameters  the  characteristics  of 
the  jammer  spectrum  must  be  specified  via  the  parameters 
of  DIGFSF.  These  parameters  are  also  members  of  NLS12. 


4.  CALLING  SEQUENCES: 

CALL  ECM  (X,  Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

5 .  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  digital  filter  routine  DIGFSF  is  used  to 
specify  the  spectral  characteristics  of  the 

j  ammer . 

b.  Flow  Chart:  Page  9-104 

c.  Cross  Reference  Table:  Page  9-221 

6.  THEORY  OF  OPERATION 


This  module  has  four  modes  of  operation: 

(1)  Spot  noise/Barrage  noise 
JPW  -  0.0 

(2)  Swept  noise  (single  sweep) 

JFMBW  /  0.0 
JPW  i  0.0 
JPEROD  =0.0 
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(3)  Swept  noise  (repeating  sweep) 

JFMBW  f  0.0 
JPW  0.0 
JPEROD  «  JPW 

(4)  Pulsed  noise 

JPW  ^0.0 
JPEROD  > JPW 

Dependeing  on  the  range  and  azimuth  positioning  of 
the  jammer  with  respect  to  the  target  either  a  self  screen¬ 
ing  or  a  standoff  jammer  can  be  simulated. 

A  block  diagram  of  this  module  is  shown  in  figure  ECM-i. 
The  output  of  the  ECM  simulation  is  summed  with  the  data  in 
arrays  X  and  Y. 


4-48 


Figure  ECM-1  BLOCK  DIAGRAM  OF  ECM 
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This  model  is  incomplete  since  the  TWT  does  not  pro¬ 
duce  saturated  noise  and  the  output  of  this  model  is  not 
properly  scaled  for  range  effects.  This  will  be  accomplished 
in  the  near  future. 
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SUBROUTINE  FGENXY 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

FGENXY  SOU-1  or  LTI-3  420 

FGENMP  SOU-1  or  LTI-3  421 

2.  PURPOSE: 

This  subroutine  is  used  to  simulate  the  radar  waveform 
generator.  In  general,  the  waveform  generator  subsystem  in 
a  radar  drives  the  power  amplifier. 

3.  INPUT  PARAMETERS 

a.  MODE  #1;  single  pulse,  internal  modulation 
Name  0/R  T  Description 


FS 

R 

F 

Simulation  sampling  rate 

F0 

R 

F 

Center  frequency  of  output  wave¬ 
form 

INORM 

R 

I 

Normalization  flag 

TI 

R 

F 

Time  increment  between  output 
waveform  samples 

NPWTX 

R 

I 

Set  --  0  for  this  mode 

FMBW 

0 

F 

Linear  FM  bandwidth 

NSUBP 

R 

I 

Set  =  0  for  this  mode 

RISTIM 

0 

F 

Rise  Time 

FALTIM 

0 

F 

Fall  Time 

TSTART 

0 

F 

Start  time  of  the  output  waveform 

PW 

R 

F 

Pulsewidth 
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Name 

o/R 

T 

Description 

VPEAK 

R 

F 

Peak  output  voltage 

N2 

R 

I 

Simulation  parameter  used  to 
specify  maximum  array  length 

Figure 

FGENXY- 1(a) 

shows  the  relationship  between 

some  of 

the 

input 

parameters  and  the  output 

waveform. 

MODE  #2;  single  pulse,  User  modulation  function 

Name  0/R  T  Description 

FS  R  F  Simulation  sampling  rate 

F0  OF  Center  frequency  of  output  wave¬ 

form 


INORM  R  I  Normalization  flag 

TI  R  F  Time  increment  between  output 

waveform  samples 

NPWTX  R  I  Number  of  points  in  user  specified 

modulation  function 


WT  R  F  Array  containing  the  user  speci¬ 

fied  modulation  function.  The 
specification  for  the  Jth  sample 
of  the  weighting  function  is  the 
following: 

WT ( 1 , J )  =  Gain 
WT(2,J)  =  Time 
WT(3,J)  =  Phase  angle 

NSUBP  R  I  Set  =  0  for  this  mode 


FSTART  0  F  Starting  frequency  at  time  = 

TSTART 


CHIRP  0  F  Linear  FM  sweep  rate 

TSTART  0  F  Start  time  of  the  output  waveform 
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Name 

O/R 

T 

Description 

VPEAK 

R 

F 

Peak  output  voltage 

N2 

R 

I 

Simulation  parameter  used  to 
specify  maximum  array  length 

NOTE: 

Either  FO  or  FSTART  &  FMBW  must  be  specified. 

Fi gure 

FGENXY-l(b) 

1  shows  the  relationship  between 

some  of 

the 

input 

parameters  and  the  output  waveform. 

MODE  #3 

;  phase  coded  waveform 

Name 

0/R 

1 

Description 

FS 

R 

F 

Simulation  sampling  rate 

F0 

0 

F 

Center  frequency  of  output 
waveform 

I  NORM 

R 

I 

Normalization  flag 

TI 

R 

F 

Time  increment  between  output 
waveform  samples 

NPWTX 

R 

I 

Set  =  0  for  this  mode 

NSUBP 

SPW 

R 

I 

Number  of  subpulses 

Subpulse  width 

PCODE 

R 

F 

Array  containing  the  phase  code 

SWTIM 

0 

F 

Switching  time  between  subpulses 

RISTIM 

0 

F 

Rise  time 

FALTIM 

0 

F 

Fall  time 

F  START 

0 

F 

Starting  frequency  at  time  =  TSTART 

CHIRP 

0 

F 

Linear  FM  sweep  rate 

TSTART 

0 

F 

Start  time  of  the  output  waveform 
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Time 


Description 


0/R  T 

VPEAK  R  F  Peak  output  waveform 

N2  R  I  Simulation  parameter  used  to 

specify  maximum  array  length 

NOTE:  Either  FO  or  FSTART  &  FMBW  must  be  specified. 

Figure  FGENXY-l(c)  shows  the  relationship  between 
some  of  the  input  parameters  and  the  output  waveform. 


4.  CALLING  SEQUENCES: 
CALL  FGENXY  (X,Y) 


Where 


X  contains  the  Output  Waveform  -  R 
Y  contains  the  Output  Waveform  -  I 


CALL  FGENMP  (X,Y) 
Where : 


X  contains  the  Output  Waveform  -  M 
Y  contains  the  Output  Waveform  -  P 


RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  number  of  waveform  samples  to  be  generated 
must  be  less  than  2**N2  where  N2  <13. 


b.  The  normalization  flag,  INORM,  is  used  to  indicate 
whether  the  output  is  to  be  interpreted  as  a  mat¬ 
ched  filter  impulse  response  or  as  a  frequency 
generator  output  waveform. 

If  INORM  =0,  2  or  4  the  output  represents  a 
frequency  generator  output. 

If  INORM  =  1  or  3  the  output  represents  the 
impulse  response  of  an  ideal  matched  filter 
and  therefore  TSTART  is  set  equal  to  zero  and 
no  trailing  zeros  are  added  to  the  output. 

This  procedure  is  required  because  the  normali¬ 
zation  procedure  used  in  the  Fourier  transform 
routine  is  different  for  waveforms  and  impulse 
responses . 

c.  For  Mode  #3  the  sign  of  NS17BP  indicates  whether 

a  polyphase  or  binary  phase  code  is  to  be  generated. 
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I, _ —  TI*(2**N2) - 

Figure  FGENXY-la  OUTPUT  WAVEFORM  ENVELOPE 


-  MODE  1 


Amplitude 


USER  SPECIFIED 
COMPLEX  ENVELOPE 


RISTIM 


FALTIM 


SWTIM 


TI*(2**N2) 


Fi2ure  FGENXY-lc  OUTPUT  WAVEFORM  ENVELOPE  -  MODE  3 


If  NSUBP  is  positive  the  phase  is  either  0°  or 
180°  with  respect  to  the  carrier.  If  NSUBP  is 
negative  the  phase  between  subpulses  is  deter¬ 
mined  by  linear  interpolation. 

d.  Operation  of  this  module  in  Mode  #3  causes  the 
phase  code  array  (PCODE)  to  be  copies  into  the 
storage  array  XMITPC  for  later  use  by  the  phase 
decoder  routine,  PHDEC.  In  this  way  the  phase 
decoder  is  automatically  slaved  to  the  phase 
code  of  the  waveform  generator. 

e.  Abort  error  codes 

1,  User  attempted  to  specify  both  a  phase  code 
and  a  complex  envelope 

2,  User  attempted  to  specify  a  complex  envelope 
with  less  than  4  points 

3,  The  time  span  of  the  waveform  to  be  generated 
is  greater  than  the  simulation  time  span, 
i.e.  TI*(2**N2). 

104,  An  error  occurred  during  execution  of  the 

weighting  function  routine,  WEITCP  or  WEITMP. 

f.  Flow  Chart:  Page  9-72 

g.  Cross  Reference  Table:  Page  9-217. 

6.  THEORY  OF  OPERATION 

This  simulation  is  used  to  simulate  the  master  frequency 
source  of  a  radar  system.  The  basic  mechanization  equation 
is  given  by  the  following  expression: 

S(t)  =  a(t)  exp  j(2rrf0t  +  0(t>) 

The  envelope  a(t)  and  the  phase  modulation  term  9(t)  are 
determined  either  directly  by  the  user  (MODE  #2)  or  indir¬ 
ectly  through  parameters  (MODE  #1  or  #2). 

This  module  can  also  handle  a  staggered  PRF  by  specify¬ 
ing  NPRIS  and  PRI(J)  in  the  system  module  (301)  namelist. 
NPRIS  is  the  number  of  pulse  repetitions  and  must  be  <11. 
PRI(J)  are  the  intervals  in  seconds  of  the  pulses.  For  ex¬ 
ample,  under  $NL301  might  be  included  NPRIS  ■  4, 

PRI  -  1.0,  0.9,  1.1,  1.0$." 
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SUBROUTINE  FILT 


MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


FILT  LTI-1  407,408 

PURPOSE : 


This  subroutine  simulates  a  continuous  filter  which 
is  defined  by  an  S-domain  polynomial  transfer  function. 


INPUT  PARAMETERS: 
Name  O/R  T 

NZ  R  I 

NP  R  I 

SF  R  F 

FZERO  0  F 

FPOLE  0  F 


Description 

Number  of  zeros  in  the  transfer 
function. 

Number  of  poles  in  the  transfer 
function. 

Scale  factor. 

Array  containing  the  zero  speci¬ 
fications.  The  s-plane  location 
for  the  Kth  zero  is  given  by  the 
following  co&fficients: 

FZERO (1,K)  =  real  component 

FZER0(2,K)  =  imaginary  component 

Array  containing  the  pole  speci¬ 
fications,  The  s-plane  location 
for  the  Kth  pole  is  given  by  the 
following  coefficients: 


FP0LE(1,K)  =  real  component 
FPOLE (2, K)  =  imaginary  component 


CALLING  SEQUENCES: 

CALL  FILT(X,Y) 

Where:  X  contains 

Y  contains 

X  contains 

Y  contains 

RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  A  purely  imaginary  pole  may  not  be  specified  within 
the  frequency  extent  of  the  input  waveform.  This 
would  cause  either  zero  or  an  extremely  small 
number  to  occur  in  the  denominator  of  the  transfer 
function. 

b.  Flow  Chart:  Page  9-55 

c.  Cross  Reference  Table:  Page  9-215 

THEORY  OF  OPERATION 


the  Input  Waveform  -  R 
the  Input  Waveform  -  I 
the  Output  Waveform  -  R 
the  Output  Waveform  -  I 


The  S-domain  transfer  function  for  a  general  filter 
can  be  represented  by  the  ratio  of  two  polynomials  as 
shown  in  the  following  expression 

(S-Z1)(S-Z2)(S-Z3).  .  .(S-ZNZ) 

H(S)  "  (S-P!)(S-P2)(S-P3).  .  .(S-Pnp)‘ 

In  general,  H(S)  is  a  complex  function  of  the  frequency 
variable  s.  The  filter  output  signal  SG(f)  can  be 
determined  using  the  following  expression 

S0(f)  -  H(f )  Si(F) 


where  Sj.(f)  is  the  input  signal  representation  in  the 
frequency  domain.  The  discrete  representation  of  these 
equations  used  in  the  module  are  the  following: 


H(J )  «  SF  . 


^K-l2!  j  FREQ  * *  CFZERe  (1 ,  K)+j  FZERO  (2 ,  K) ]  ^ 
/7K*1P{  j  FREQ  3  )  “  C FP6LE  ( 1 ,  K )+j  FP0LE  ( 2 ,  K) ]  } 


X(J)+j  Y(J)  =  H(J)  *  [x(J)+j  Y(J)] 


where:  FREQ(J)  is  the  frequency  associated  with  the 

Jth  sample 
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SUBROUTINE  HET 


1.  MODULE  IDENTIFICATION: 


Name 


Classification  Code  Reference  Number 
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Figure  HET-1  BLOCK  DIAGRAM  OF  SSB  MODULATOR 
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SUBROUTINE  HWDET 


1 .  MODULE  IDENTIFICATION : 


Name 

Classification  Code 

Reference  Number 

HWDET 

NL-1 

414,415 

FWDET 

NL-1 

416,417 

SQDET 

NL-1 

418,419 

HLIM 

NL-1 

438,439 

IHWDET 

NL-1 

445,446 

IFWDET 

NL-1 

447,448 

I  SQDET 

NL-1 

449,450 

IHLIM 

NL-1 

442,443 

2,  PURPOSE: 

This  subroutine  simulates  a  signal  detector . 
Depending  on  the  entry  point,  a  half-wave  detector, 
a  full-wave  detector,  a  square-law  detector,  or  a 
hard  limiter  is  simulated,, 

3 „  INPUT  PARAMETERS: 

NONE 

4„  CALLING  SEQUENCES  AND  THEORY  OF  OPERATION 
a.  Half-Wave  Detector 

CALL  HWDET (XIN,XOUT) 

Where:  XIN  contains  the  Input  Waveform 

XOUT  contains  the  Output  Waveform 
XOUT(J)  =  XIN(J)  ;  if  XIN(J)  >  0,0 
=  0,0  ;  otherwise 
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b.  Full-Wave  Detector 
CALL  FWDET (XIN , XOUT ) 

Where:  XIN  contains  the  Input  Waveform 

XOUT  contains  the  Output  Waveform 

XOUT(J)  =  |  XIN ( J ) | 

c.  Square-Law  Detector 
CALL  SQDET (XIN, XOUT) 

Where:  XIN  contains  the  Input  Waveform 

XOUT  contains  the  Output  Waveform 
XOUT(J)  =  XIN  (J  )*XIN  (J  ) 

d.  Hard  Limiter 

CALL  HLIM (XIN, XOUT) 

Where:  XIN  contains  the  Input  Waveform 

XOUT  contains  the  Output  Waveform 
XOUT(J)  =  1.0  ;  if  XIN(J)  >  0.0 
=  -1.0;  if  XIN(J)  <  0.0 

e.  Digital  Half-Wave  Detector 
CALL  IHWDET (IN, IOUT ) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 
IOUT(J)  -  IOUT(J)  ;  if  IOUT(J)  >  0 
■  0  ;  otherwise 
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f.  Digital  Full-Wave  Detector 
CALL  IFWDET(IN,IOUT) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 
IOUT(J)  =  IABS(IN(J>) 

g.  Digital  Square-Law  Detector 
CALL  ISQDET(IN, IOUT) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 
IOUT(J)  -  IN(J)*IN(J) 

h.  Digital  Hard  Limiter 
CALL  IHLIM ( IN , IOUT ) 

Where:  IN  contains  the  Input  Waveform 

IOUT  contains  the  Output  Waveform 
IOUT(J)  -  1  ;  IN ( J )  >  0 

-  0  ;  IN(J )  **  0 

-  1  ;  IN  ( J  )  <  0 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-137 

b.  Cross  Reference  Table:  Page  9-225 


4-67 


y.tA 

SUL-t.i’lil  1M-  tiwl'l  Tl  xis.  xt  tu  ) 

..Ki-NsKN  X1NI  ll.XUUTIll.INin.IuUTIll 

yjt>s 

UAIrt  «L,/ 

>5^-5 

NS  -  Inl.uLlXlN(Nl^3l> 

3">o^ 

LU  oU  1-1 t  MS 

XltUI  111  -  X1N(  1  1 

S  >c-c> 

>■  IK  XLtil  til  -Li.  °*°  1  XOUTtll  =  ‘'•u 

J*>  o  / 

1,1.  In  K.u 

l 

i  » .  ,v 

rNIKY  Kl-t  1  <  X 1  N  *  XlUTI 

>  .»  / 

NS  -  ibt.L'L  (X1N(N1w-4)1 

yj'/i 

bL  **u  A- 1 

iji. 

nl  fill  s  Aoi  (X  IN(  i  1  1 

1  L  \  U  1 

v/f- 

L 

...UY  sioiH  xin,  xlui) 

>*>  Yv> 

Ns  -  Ita'lllihl1*!'’-'  >1 

' 

_<  » /  # 

I'U  *■  1  1  i  -  1  t  f*«> 

•>  r  j 

ol,  ,.,,,itii  -  Al'.(ii*xl:Hll 

••> 

i,i  n,  Ijo 

.>*  .  u 

U 

.  Ax.* 

tMKI  I.L  1 1  IAIN.  XH'1  1 

/  V  * 

U.  -  1  Kl'tll  (  A  JNINl  vl  )  ) 

, 

,n.  o.i  1-1. NS 

ll  (  X  1  N 1  i  •  )  S 1  ,  1  /  »  1’ . 

,  >U  )  (  1  I--1.U 

s  ., 

l.l  U  Ui 

J.*V,  l 

Al-l.l  (  1  1  =  1  •' 

r  (  u 1 1 nu 

>  , ,  - 

i. 

,  .  |  |  -  XI  M.’ll'i  '  1 

>•  •  1 

XI  1. 1  1  Nl'i  »  1  -  XIMNK-  » 

...ll  INKt.l  =  XI  MM'Y  1 

iil.il  N 

'  •' K*  •* 

l 

t  V  It  N  lit  l  4  |.M.  1 1'*n  » 


SUBROUTINE  INGTOR 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

INGTOR  LTI-4  409,  410 

INGNCL  LTI-4  None 

2.  PURPOSE: 

This  subroutine  is  used  to  simulate  a  digital  inte¬ 
grator  . 

3.  INPUT  PARAMETERS : 

Name  0/R  T  Description 

FBCK  R  F  Feedback  coefficient 

4.  CALLING  SEQUENCES: 

CALL  INGTOR  (DIN,  DOUT) 

Where:  DIN  contains  the  Input  Waveform 

DOUT  contains  the  Output  Waveform 

The  storage  register  (Cl)  is  cleared  before 
execution  begins. 

CALL  INGNCL  (DIN,  DOUT) 

Where:  DIN  contains  the  Input  Waveform 

DOUT  contains  the  Output  Waveform 

The  storage  register  (Cl)  is  not  cleared 
before  execution  begins. 

5.  RESTRICTIONS,  REQUIREMENT.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-97 

b.  Cross  Reference  Table:  Page  9-220 
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6. 


THEORY  OF  OPERATION 


The  block  diagram  of  the  digital  integrator  simulated 
by  this  module  is  shown  in  Figure  INGTOR-1.  The  Z-plane 
transfer  function  is  given  by  the  following  expression: 

T(Z)  "  Z-FBCK 
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SUBROUTINE  IONOS 


I.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


IONOS 


LTI-1 


511 


2.  PURPOSE: 


This  subroutine  is  used  to  simulate  the  effect  of 
the  ionosphere  on  signals  propagating  through  it. 


3.  INPUT  PARAMETERS: 
Name  0/R  T 

RFF0  R  F 


SEDENS  R  F 


Description 

Center  frequency  of  the  electro¬ 
magnetic  wave  traversing  the 
ionosphere 

The  integrated  electron  density 
along  the  propagation  path 


4.  CALLING  SEQUENCES: 

CALL  IONOS  (X,Y) 

the  Input  Waveform  -  R 
the  Input  Waveform  -  I 
the  Output  Waveform  -  R 
the  Output  Waveform  -  I 


Where:  X  contains 

Y  contains 
X  contains 

Y  contains 


5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  effect  of  the  earth’s  magnetic  field  is  not 
included  in  this  model. 
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bo  The  effect  of  electron  collisions  with  ions  and 
neutral  particles  is  not  included  in  this  models 
The  major  effect  resulting  from  electron  collisions 
is  a  damping  term  which  results  in  absorption  of 
energy  from  the  electromagnetic  wave  transversing 
the  ionized  medium0  For  frequencies  above  100  MHz 
electron  collisions  have  negligible  effect  on  the 
phase  behavior  of  the  ionosphere  as  a  function  of 
frequency, 

c  .  Reference : 

Berkowitz,  R3  S,  ed.,  Modern  Radar  Analysis. 
Evaluation,  and  System  Design,  John  Wiley  &  Sons, 
Inc c ,  New  York,  1965,  pp  349-354,  364-370, , 

d  Flow  Chart:  Page  9-117 

ec  Cross  Reference  Table:  Page  9-223 

6  THEORY  OF  OPERATION 


The  index  of  refraction  of  an  ionized  medium  is 
given  by  the  following  expression: 


n 


Ne  e~ 

_c2 

mml: 


(1) 


where : 


e 


f 


is  the  electron  density  (electrons /crrr  ) 

is  the  electron  charge  (48  x  10“^®  esu) 

is  the  frequency  of  this  incident  electro¬ 
magnetic  energy 


Collisions  between  electrons  and  ions  or 
neutral  particles  are  neglected. 


The  phase  velocity  of  the  propagating  wave,  Vp,  is 
given  by  the  following  expression: 

Vp  =  c/n  (2) 

where:  c  is  the  speed  of  light  in  free  space 
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Therefore, combining  equations  1  and  2  the  following 
expression  is  obtained: 


(3) 


The  phase  velocity  will  approach  infinity  as  the 
denominator  of  equation  3  decreases  to  zero.  When 
this  condition  arises,  further  wave  propagation  is 
impossible,.  The  frequency  for  which  this  occurs  is 
called  the  critical  frequency  and  is  given  by  the 
following  expression: 


9  np  2 

fr  = - -  =  8.0  x  10'  NP  (Hz) 


m  7t 


(4) 


Equation  3  now  becomes  the  following: 
C 


VP  * 


\h  -  fc2/f2 


(5) 


The  differential  phase  shift  d0  that  a  wave  will 
encounter  in  traversing  an  element  of  path  length  ds 
is  given  by  the  following  expression: 


(6) 


where:  \  is  the  wavelength  of  the  radiated  electro¬ 

magnetic  energy 

Since  Vp/f  the  following  expression  is  obtained 
by  substitution  of  equation  5  into  equation  6: 


d0 


2mf  / 

C  \/  1  " 


f  c  2  /  f  2 


ds 


(7) 


This  expression  can  be  further  simplified  by  adding 
the  term  fc4  to  the  terms  inside  the  square  root 

4f2 
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radicals  For  most  radar  applications  f c  «  f  so  this 


will  introduce  a  very  small 
becomes : 


error „  Equation  7  now 


ds 


ds  (8) 


In  order  to  determine  the  phase  shift  over  a  propagation 
path  from  si  to  S2  the  following  expression  will  be 
used 


2 

substituting  for  f  the  following  is  obtained: 


0  = 


ds 


Cf 


8.0  x  10 


SJ 


Ne  ds 


(9) 


The  first  term  of  equation  9  is  the  linear  phase  shift 
due  to  propagation  over  a  path  of  length  S2  -  s[  and 
wi 1 l  be  omitted  since  it  represents  a  time  delay  onlv 

Therefore,  ihe  equation  used  to  calculate  the 
disperive  effect  of  the  ionosphere  is  given  by  the 
following  expression: 


4-7  7 


0(f)  =  “2  TT 

where:  SEDENS 

f 

FEXT 

F 


—  ' 0  -  \  SEDENS 

3o0  x  1010  f 


integrated  electron 
density 


FEXT  FEXT 

”2  AF  k  '  2  AF 

frequency  extent  of  input  waveform 
representations 

independent  variable  spacing  of  input 
waveform  representations „ 


s2 

Ne  ds  = 

S1 

RFF0  +  k  AF, 
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SUBROUTINE  LAMPCP 


MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

LAMPCP  LTI-1  or  LTI-4  458 

LAMPRE  LTI-1  or  LTI-4  456,  457 

PURPOSE: 

This  subroutine  is  used  to  simulate  a  linear 
amplifier. 

INPUT  PARAMETERS: 

Name  O/R  T  Description 

GAIN  R  F  Amplifier  gain 

CALLING  SEQUENCES: 

CALL  LAMPCP  (XIN,  YIN,  XOUT,  YOUT) 

Where:  XIN  contains  the  Input  Waveform  -  R 

YIN  contains  the  Input  Waveform  -  I 
XOUT  contains  the  Output  Waveform  -  R 
YOUT  contains  the  Output  Waveform  -  I 
CALL  LAMPRE  (XIN,  XOUT) 

Where:  XIN  contains  the  Input  Waveform 

XOUT  contains  the  Output  Waveform 
RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-68 

b.  Cross  Reference  Table:  Page  9-217 
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6.  THEORY  OF  OPERATION 


The  relationship  between  the  input  and  output  is 
given  by  the  following  expressions: 


XOUT(J) 

YOUT(J) 


GAIN  *  XIN(J)  ;  LAMPCP  &  LAMPRE 
GAIN  *  YIN ( J )  ;  LAMPCP 
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SUBROUTINE  MTIFLT 


I.  MODULE  IDENTIFICATION: 
Name 


MTIFLT  (MODEDF=I) 
MTINCL  (MODEDF-1) 
MTIIFT  (MODEDF=2 ) 
MTIINC  (MODEDF=2 ) 
2.  PURPOSE: 


Classification  Code 
LTI-4 
LTI-4 
NL-2 
NL-2 


Reference  Number 
430,431 
432,433 
430,431 
432,433 


This  subroutine  is  used  to  simulate  a  double  delay 
MTI  processor.  Either  floating  point  or  integer  arithmetic 
is  selectable. 


3.  INPUT  PARAMETERS 


a.  MODEDF  =  1 


Name 

0/R 

T 

Description 

FF0 

0 

F 

Feed  -  forward  coefficient  -  0  dela 
(no  delay) 

FFl 

0 

F 

Feed  -  forward  coefficient  -  1  dela 
(single  delay) 

FBI 

0 

F 

Feedback  coefficient  -  1  delay 

FB2  0 

b .  M0DEDF=2 

F 

Feedback  coefficient  -  2  delay 
(double  delay) 

Name 

0/R 

T 

Description 

IFF0N 

0 

I 

Numerator  of  Feed  -  forward  coef  - 
0  delay 

IFF0D 

R 

I 

Denominator  of  feed  -  forward  coef- 
0  delay 

IFF1N 

0 

I 

Numerator  of  feed  -  forward  coef  - 
1  delay 

Name 


0/R  T 
IFF1D  R  I 

IFB1N  0  I 

IFB1D  R  I 

IFB2N  0  I 

IFB2D  R  I 


Description 

Denominator  of  feed  -  forward  coef  - 
1  delay 

Numerator  of  feedback  coef  -  1  delay 

Denominator  of  feedback  coef  - 

1  delay 

Numerator  of  feedback  coef  -  2  delay 

Denominator  of  feedback  coef.  - 

2  delay 


NBITS  R  I  Number  of  bits  to  be  used  in 

storing  one  range  sample  in  the 
arrays  (ISR1  and  ISR2)  which 
represents  the  digital  delay 
lines  (includes  sign  bit). 


CALLING  SEQUENCES: 


Floating  point  arithmetic 


CALL  MTIFLT  (DIN,  DOUT) 


Where:  DIN  contains  the  Input  Waveform 


DOUT  contains  the  Output  Waveform 

The  range  bit  storage  arrays  (SRI  and  SR2) 
are  cleared  prior  to  execution. 

The  entry  point  is  used  only  if  M0DEDF=1. 

Floating  point  arithmetic 

CALL  MTINCL  (DIN,  DOUT) 

Where:  DIN  contains  the  Input  Waveform 

DOUT  contains  the  Output  Waveform 

The  range  bin  storage  arrays  (SRI  and  SR2) 
are  not  cleared  prior  to  execution  except 
for  the  first  execution  of  this  subroutine. 
This  entry  point  is  used  only  if  MODEDF-l. 


Integer  arithmetic 


CALL  MTII FT  ( LN ,  I OUT) 

Where:  LN  contains 

IOUT  contains 

The  range  bin 
cleared  prior 

This  entry 


the  Input  Waveform 
the  Output  Waveform 
arrays  (1SR1  ami  I S K 2 )  are 


to  execut ion . 
point  is  used  only  if  M0DEDF=2 . 


Integer  Arithmetic 


CALL  MTIINC  (IN,  OUT) 

Where:  LN  contains  the  input  Waveform 

IOUT  contains  the  Output  Waveform 

The  range  bin  storage  arrays  (ISRl  and  ISR2) 
are  not  cleared  prior  to  execution  except  for 
the  first  execution  of  this  subroutine. 


This  entry  point  is  used  only  if  M0DEDF=2 . 

5.  RESTRICTIONS,  REQUIREMENTS ,  MISCELLANEOUS  DATA 

a.  The  maximum  number  of  range  bins  simulated  is  2048. 
If  more  than  2048  range  samples  are  contained  in  the 
input  array,  only  the  first  2048  are  processed. 

b.  The  parameter  MODEDF  is  used  by  the  simulation  con¬ 
troller  to  determine  whether  the  floating  point  or 
the  fixed  point  module  entry  points  are  to  be  used. 

c.  The  entry  points  which  use  integer  arithmetic  are 
considered  nonlinear  because  of  truncation  in  arith¬ 
metic  and  saturation  in  the  delay  line. 

d.  Flow  Chart:  Page  9-201 

e.  Cross  Reference  Table:  9-234 


4-85 


THEORY  OF  OPERATION 


The  block  diagram  of  the  two-delay  MTI  processor  digital 
filters  simulated  by  this  module  are  shown  in  Figure  MTIFLT-1 
and  MTIFLT-2.  The  Z-plane  transfer  function  for  each  range 
bin  (floating  point  arithmetic)  is  given  by  the  following 
expression : 

2  FF1  1 
2  -I-  TT0  Z  + 

T(2)  * FF0  z2  :~FB1  TT^r 


The  Z-plane  transfer  function  for  each  range  bin  (integer 
arithmetic)  is  approximated  by  the  following  expression: 

2  IFF1N  .  IFF0D  _  IFF0D 

_  IFF0N  zL  +  IFF0N  •  IFF1D  IFF0N 

T  Z  IFF0D  2  ifbIN  IFB2N 

Z  "  IFBlD  Z  *  IFB2D 


The  delay  represented  by  the  Z  operator  is  determined  by  the 
radar  pulse  repetition  interval. 


Figure  MTIFLT-1  BLOCK  DIAGRAM  OF  MTIFLT/MTINCL 

(Floating-point  Arithmetic) 


IFF0N 

IFF0D 


Figure  MTIFLT-2  BLOCK  DIAGRAM  OF  MTIFLT/MTINCL  (Integer  Arithmetic) 
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SUBROUTINE  NONLIN 


I*  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


NONLIN 
20  PURPOSE : 


NL-1 


401,  402 


This  subroutine  is  used  to  simulate  a  zero  memory 
nonlinear  device „ 


3o  INPUT  PARAMETERS: 
Name  O/R  T 

NPTS  R  I 

TFN  R  F 


Description 

Number  of  points  used  to  specify 
the  transfer  function 

Array  containing  the  specification 
of  the  transfer  function,,  The 
coefficients  for  the  Jth  point 
are  the  following: 


TFN(1,J)  =  Input  voltage 
TFN(2,J)  =  Output  voltage 


4 .  CALLING  SEQUENCES : 

CALL  NONLIN  (A,$mmmm) 

Where:  A  contains  the  Input  Waveform 


A  contains  the  Output  Waveform 


$mmmm  is  the  statement  number  in  the 

calling  program  to  which  control 
is  returned  if  a  discrepancy  is 
detected  in  the  input  data. 
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5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  independent  variable  in  the  transfer  function 
specification  array  must  be  monotonically  increasing. 

b.  The  input  data  must  be  within  the  limits  of  the 
specified  transfer  function. 

c.  Flow  Chart:  Page  9-135 

d.  Cross  Reference  Table:  Page  9-225 

6.  THEORY  OF  OPERATION 


The  input  vs  output  characteristics  of  the  device  to  be 
simulated  are  specified  via  the  array  TFN.  A  table 
lookup  and  linear  interpolation  scheme  are  used  to  com¬ 
pute  the  output  value  for  each  sample  of  the  input 
waveform.  An  abnormal  termination  will  occur  if  an 
input  sample  is  outside  the  range  of  input  values  speci¬ 
fied  in  the  array  TFN. 

An  example  of  the  input  data  for  this  module  is  the 
nonlinear  transfer  function  shown  in  Figure  NONL IN- 1 . 

The  associated  input  data  statement  is  as  follows: 

$NL401  TFN  =  -2.0,  -2.0,  -1.0,  -2.0,  1.0,  2.0,  2.0, 

2.0,  NTFN  =  4$ 
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OUTPUT 


Figure  NONLIN-1  EXAMPLE  OF  NON-LINEAR  TRANSFER  FUNCTION 
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SUBROUTINE  PHDEC 


1 .  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


PHDEC  LTI-2  508,  509 

2.  PURPOSE: 

This  subroutine  is  used  to  simulate  an  analog  phase 
decoder . 

3 .  INPUT  PARAMETERS 


a.  Automatic  mode;  INPTF  =  0 


Name  O/R  T 

SIMF0  R  F 

INPTF  R  T 

SPW  R  F 

XMITPC  R  F 

NSUBP  R  I 


Description 


Center  frequency  of  waveform 
previously  generated  by  FGENXY/ 
FGENMP.  This  parameter  is  used 
to  specify  the  center  frequency 
of  the  tapped  delay  line. 

Set  =  0  for  this  mode 

Subpulsewidth  of  waveform  pre¬ 
viously  generated  by  FGENXY/ 
FGENMP.  This  parameter  is  used 
to  specify  the  delay  line  tap 
spacing  (TAPSPC). 

Array  containing  the  phase  code 
of  the  waveform  generated  by 
FGENXY /FGENMP. 

Number  of  subpulses  in  the  wave¬ 
form  generated  by  FGENXY/FGENMP . 
This  parameter  specifies  the  num¬ 
ber  of  delay  line  taps. 


NOTE:  No  user  supplied  inputs  required  for  this  mode. 

All  required  parameters  are  initialized  by 
FGENXY/FGENMP. 
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b.  User  Supplied  data  mode;  INPTF  =  1 


Name  O/R  T 

F0DEC  R  F 

TAPSPC  0  F 

INPTF  R  I 


Description 

Center  frequency  of  the  tapped 
delay  line. 

Spacing  between  delay  line  taps 
This  parameter  is  set  equal  to 
SPW  if  it  is  not  specified  by 
the  user. 

Set  =  1  for  this  mode 


XMITPC  R  F  Array  containing  the  phase  code 

of  the  waveform  generated  by 
FGENXY/FGENMP.  (Initialized  by 
FGCENXY/FGENMP) 

NSUBP  R  I  Number  of  subpulses  in  the  wave 

form  generated  by  FGENXY/FGENMP 
This  parameter  specifies  the 
number  of  delay  line  taps. 
(Initialized  by  FGENXY /FGENM) 

4.  CALLING  SEQUENCES: 

CALL  PHD EC  (X,Y) 


Where:  X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 


5 •  RESTRICTIONS ,  RECOMMENDATIONS,  MISCELLANEOUS  DATA 

a.  If  INPTF=1  and  TAPSPC  is  specified,  then  TAPSPC 
should  be  an  integral  multiple  of  the  period  of 
the  waveform  center  frequency,  i.e.  TAPSPC  = 
N*1.O/SIMF0  where  N  is  an  integer.  Otherwise, 
the  taps  will  not  have  the  desired  phase  shifters 
as  specified  in  the  array  XMITPC. 

b.  Flow  Chart:  Page  9-78 

c.  Cross  Reference  Table:  Page  9-218 
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h.  THEORY  OF  OPERATION 

This  module  generates  the  transfer  function  of  a  tapped 
analog  delay  line.  The  main  use  of  this  module  is  to  simu¬ 
late  a  surface  acoustic  wave  device.  A  block  diagram  of  this 
device  represented  by  this  module  is  shown  in  Figure  PHDEC-1. 
This  module  is  structured  to  represent  only  the  taps  (energy 
pickoffs)  of  the  delay  line.  The  wave  launches  must  be 
specified  separately  in  terms  of  its  transfer  function. 

This  can  be  done  via  either  the  FILT  or  the  WEIT  modules. 

A  pictorial  diagram  of  a  surface  acoustic  wave  device  and 
its  representation  are  shown  in  Figures  PHENC-2(a)  and  -2(b). 
The  impulse  response  of  the  device  represented  by  this  module 
is  given  by  the  following  expression: 

K-NSUBP 

h(t)-2j  gK  S(t-Ktj)  eJ  " 

K-l 

where:  gK  =  taP  gai-n 

t^  =  tap  spacing 
0K  =  tap  phase  shift 

The  transfer  function  computed  by  this  module  is  given  by 
the  following  expression: 

rx  NSUBP 

n  ,  /  .  V1  j  0k  -j2mft 

H  (f )  =  |  [  h  ( t )]  =#  2Lt  \  8(t-ktd)e  e  dt 

-x  K  1 

Interchanging  the  order  of  integration  and  summation 
NSUBP  . 

h ( f )  =  X  §keJ 
1 

NSUBP 

Ej  0k  -  j  2  rrf  k  t . 

«ke  e 

1 

This  version  of  the  phase  decoder  module  sets  all  g^  =  1.0. 
if  a  phase  decoder  with  gain  weighting  or  errors  is  desired 
then  the  DFT  routine  should  be  used  to  compute  the  transfer 
function . 
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Figure  PHDEC-1  BLOCK  DIAGRAM  OF  A  TAPPED  DELAY  LINE  PHASE  DECODER 
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SUBROUTINE  PHENC 


1. 

MODULE  IDENTIFICATION: 

Name 

Classification  Code 

Reference  Number 

PHENC /FGGNXY 

SOU-1  or  LTI-3 

506 

PHENC /FGENMP 

SOU-1  or  LTI-3 

507 

2  .  PURPOSE: 

This  subroutine  generates  a  phase  code  which  is  subse¬ 
quently  used  by  FGENXY/FGENMP  to  generate  a  binary  phase 
coded  waveform.  Either  the  module  FGENXY  or  FGENMP  is 
automa tica lly  scheduled  for  execution. 

3.  INPUT  PARAMETERS 


a.  MODEPH  =  1;  Barker  Phase  Code 


Name 

MODEPH 

NSUBP 


0/R 

R 

R 


I 

I 

I 


Description 

Set  =  1  for  this  mode 

Number  of  subpulses  to  be  generated. 
The  allowable  values  of  this  para¬ 
meter  and  the  corresponding  Barker 
codes  are  as  follows: 

2  +  - 

3  +  +  - 

4  +  +  -  4 

5  4  +  4-  -  4- 

7  4-  4-  4-  -  -  4-  - 

11  4-4-4-  -  --4--  -  4-  - 

13  +  +  +  +  +--++-  +  -  + 


b.  MODEPH  =  2; 

Pseudo  Random  Phase  Code 

Name 

0/R  1 

Description 

MODEPH 

R  I 

Set  =  2  for  this  mode 

NSUBP 

R  I 

Number  of  subpulses, 
value  -  300  for  this 
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UNCLASSIFIED 


Name 


0/R  T 


IPY  R  I 


NSR  R  I 


Description 

Array  representing  the  shift 
register  used  to  generate  the 
pseudo  random  sequence.  The 
storage  cells  of  this  shift 
register  contains  either  1  or  0. 
The  user  can  specify  the  initial 
condition  of  this  register  if 
desired. 

Number  of  stages  in  the  shift 
register. 


c.  MODEPH  =  3;  User  Specified  Phase  Code 


Name  0/R  T  Description 

MODEPH  R  I  Set  =  3  for  this  mode 

NSUBP  R  I  Number  of  subpulses.  For  this 

mode  this  parameter  is  limited 
to  35. 

CODE  R  B  Word  containing  the  phase  code. 

Each  bit  of  this  36  bit  word 
represents  the  phase  of  one  sub¬ 
pulse.  The  phase  for  the  first 
subpulse  is  determined  by  the 
LSB;  i.e.  right  most  bit. 


CALLING  SEQUENCES: 
CALL  PHENC  ($mmmm) 


Where:  mmmm  is  the  statement  number  to  which 

control  will  be  transferred  if  an  error 
is  detected  in  the  input  data. 


RESTRICTIONS.  RECOMMENDATIONS.  MISCELLANEOUS  DATA 

a.  Reference  for  Barker  and  Pseudo  random  phase  codes 

Skolnik,  M.  I.:  RADAR  HANDBOOK.  McGraw-Hill  Book 
Company,  New  York,  1970,  pp  20-19  and  20-20. 

b.  Flow  Chart:  Page  9-167 

c.  Cross  Reference  Table:  Page  9-229 
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THEORY  OF  OPERATION 


This  module  generates  a  phase  code  determined  by  the 
input  parameters  and  loads  the  array  PCODE  with  the  phase 
code.  The  operations  performed  for  MODEPH  *  1  or  3  consist 
of  connecting  bits  of  the  word  CODE  into  phase  angles.  A 
binary  "1"  is  converted  into  180  degrees  and  a  binary  "0"  is 
converted  into  0  degrees. 

Figure  PHDEC-1  is  a  block  diagram  of  the  shift  regis¬ 
ters  used  to  generate  the  pseudo  random  phase  codes.  The 
parameter  NSR  determines  the  selection  of  shift  register 
taps  (Reference  Skolnik,  p  20-20,  table  6).  The  output 
from  the  delay  line  is  converted  to  a  phase  code  according 
to  the  procedure  outlined  in  the  preceeding  paragraph. 
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SUBROUTINE  RDIGFL 


1. 

MODULE 

IDENTIFICATION: 

Name 

Classification  Code 

Reference  Number 

RDIGFL 

LTI-4 

405 

RDFNCL 

LTI-4 

406 

2.  PURPOSE; 

This  subroutine  is  used  to  simulate  a  double  delay 
digital  filter  and  is  capable  of  synthesizing  a  filter  trans¬ 
fer  function  having  two  poles  and  two  zeros. 


3.  INPUT  PARAMETERS ; 


Name 

0/R 

T 

Description 

FF0 

0 

F 

Feed  -  forward  coefficient  -  0  delay 

FF1 

0 

F 

Feed  -  forward  coefficient  -  1  delay 

FBI 

0 

F 

Feedback  coefficient  -  1  delay 

FB2 

0 

F 

Feedback  coefficient  -  2  delay 

4. 


CALLING  SEQUENCES: 


CALL  RDIGFL  (X,Y) 


Where : 


X  contains 

Y  contains 
X  contains 

Y  contains 


the  Input  Waveform  -  R 
the  Input  Waveform  -  I 
the  Output  Waveform  -  R 
the  Output  Waveform  -  I 


The  storage  registers  (RR1,  RR2,  RI1,  and 
RI2)  are  cleared  before  execution  begins. 

CALL  RDFNCL  (X,Y) 


Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  output  Waveform  -  I 
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The  storage  registers  are  not  cleared 
before  execution  begins. 


5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-100 

b.  Cross  Reference  Table:  Page  9-221 

6.  THEORY  OF  OPERATION 


The  block  diagram  of  the  two  delay  digital  filter  simu¬ 
lated  by  this  module  is  shown  in  Figure  RDIGFL-1.  The  Z- 
plane  transfer  function  is  given  by  the  following  expression 

Z2+  Z  +  1/FF0 
T(Z)  =  FF0  FFQ _  . 

Z2  -  FBI  Z  -  FB2 


} 
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SUBROUTINE  RNDARY 


MODULE 

IDENTIFICATION: 

Name 

Classification  Code 

Reference  Number 

RNDARY 

SOU-1 

214,  215 

ADDRND 

SOU- IS 

235,  236 

ADRNDC 

S0U-1S 

237 

2.  PURPOSE: 

This  subroutine  loads  an  array  with  random  numbers 
generated  by  the  randoin  number  generator  function  sub¬ 
routine  RRAND. 

3.  INPUT  PARAMETERS: 


1 

Description 

NPTS 

A 

I 

Number  of  random  distribution 
samples  to  be  loaded  into  the 
array  (RNDARY  only). 

NTYPE 

R 

I 

Code  indicating  the  type  of 
random  distribution  to  be  used 
(RNDARY  and  ADDRND  only). 

TI 

R 

F 

Simulation  sampling  increment 
(RNDARY  only). 

4.  CALLING  SEQUENCES: 

CALL  RNDARY(RND) 

Where:  RND  contains  the  Output  Waveform 

CALL  ADDRND(RND) 

Where:  RND  contains  the  Input  Waveform 

RND  contains  the  Output  Waveform 
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CALL  ADRNDC  (RND,  RNDY) 


Whe  re : 


RND  contains 
RNDY  contains 
RND  contains 
RNDY  contains 


the  Input  Waveform  -  R 
the  Input  Waveform  -  I 
the  Output  Waveform  -  R 
the  Output  Waveform  -  I 


This  module  is  used  only  to  add 
Gaussian  samples  to  the  input  waveform 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  User  supplied  function  subprogram  RRAND  is  called 

by  the  subprogram  and  must  be  initialized  via  NL101 

or  NL102  before  this  subroutine  is  executed. 

b.  The  random  distribution  code  (NTYPE)  is  as  follows: 

(1)  Floating  point  uniform  distribution  with 
specified  mean  and  extent. 

(2)  Floating  point  Rayleigh  distribution  with 
specified  standard  deviation. 

(3)  Floating  point  Gaussian  distribution  with 
specified  mean  and  standard  deviation. 

(4)  Positive  integer  uniform  distribution  from 
0  to  2**36. 

(5)  Not  used  by  this  subprogram. 

(6)  Swerling  Target  Models  #1  &  #2  with  specified 
mean  cross  section. 

(7)  Swerling  Target  Models  #3  &  #4  with  specified 
mean  cross  section. 

(8)  Sine  distribution. 

c.  Flow  Chart:  Page  9-150 

d.  Cross  Reference  Table:  Page  9-227 
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6 „  THEORY  OF  OPERATION 


The  mechanization  equations  for  each  entry  point 
are  the  following: 

RNDARY:  RND(J)  =  RRAND  (NTYPE)  ;  KJ  <  NPTS 

ADDRND:  RND(J)  =  RND(J)  +  RRAND  (NTYPE);  the  range  of 
J  is  determined  by  the  input  array „ 

ADRNDC:  RND(J)  =  RND(J)  +  RRAND(3) 

RNDY(J)  =  RNDY(J)  +  DUM 

where:  RRAND(3)  contains  the  cosine  gaussian 

distribution  sample 

DUM  contains  the  sine  gaussian 

distribution  sample 
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SUBROUTINE  SHIFT 


1.  MODULE  IDENTIFICATION: 


Name 

Classification  Code 

Reference  Number 

SHIFT 

LTI-1 

224,225 

SHIFTS 

LTI-1S 

226,227 

RSHIFT 

LTI-1 

229,230 

RSHFTS 

LTI-1S 

231,232 

2.  PURPOSE: 

This  subroutine  is  used  to  delay  a  waveform  in  time 
and/or  introduce  a  constant  phase  shift. 

3.  INPUT  PARAMETERS: 


Name 

0/R 

T 

Description 

T0 

R 

F 

Time  delay 

THT 

R 

F 

Phase  shift 

TJIT 

0 

F 

Time  delay  jitter 

SIMF0  R  F 

CALLING  SEQUENCES: 

Center  frequency  to  be  used 
computing  time  delay  effect 

CALL  SHIFT  (X,Y,A,B) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
A  contains  the  Output  Waveform  -  R 
B  contains  the  Output  Waveform  -  I 
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The  basic  mechanization  equation  for  this  entry  point  is 
the  following: 


A( J)+j B( J)™  [x(J)+jY(J)]*eJ[THT-2”  (FSTRT+n*J)T«] 

Where:  FI  =  frequency  increment 

FSTRT  -  frequency  corresponding  to  the 

first  element  of  the  input  array 


CALL  SHIFTS  (X,Y,A,B) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

A  contains  the  Output  Waveform  -  R 

B  contains  the  Output  Waveform  -  I 


The  basic  mechanization  equation  for  this  entry  point  is 
the  following: 

j  [THT-2tt(FSTRT  +FI*J)T0] 

A(J)+jB(J)=A(J)+jB(J)  +  lx(J)+jY(J)J  *e 


CALL  RSHIFT  (X,Y,A,B) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
A  contains  the  Output  Waveform  -  R 
B  contains  the  Output  Waveform  -  I 

The  basic  mechanization  equation  for  this  entry  point  is 
the  following: 


A(J)  ♦  3B(J)  -  [x(J)+jY(J)J  *e3  [ THT-2ir( SIMF0+FI*J)  (T04p •  TJIT)] 

where  p  is  a  sample  of  a  random  process  having  a  uniform 
distribution  from  -1/2  to  1/2. 


CALL  RSHFTS(X , Y ,A , B) 


Where : 


X  contains  the  Input  Waveform  -  R 
Y  contains  the  Input  Waveform  -  I 
A  contains  the  Output  Waveform  -  R 
B  contains  the  Output  Waveform  -  I 


The  basic  mechanization  equation  for  this  entry  point  is 
the  following: 


A(J)+jB(J)=A(J)+jB(J)  +  |x(J)+jY(J)]*e 


j [THT-2ff(SIMF0+FI*J) *(T0+p.TJIT3 


5.  RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-63 

b.  Cross  Reference  Table:  Page  9-216 

6.  THEORY  OF  OPERATION 


s(t)  : 


If  the  Fourier  transform  pair  is  given  for  the  signal 


s(t) 


S(f) 


then  s(t-tQ) 


S(f)  e'j2,rfto 


This  expression  was  derived  in  reference  1. 

In  addition  the  introduction  of  a  constant  phase  shift 
in  the  time  waveform  results  in  the  same  constant  phase 
shift  in  the  frequency  spectrum,  e.g. 


s(t)e- 


S(f)  e“ 


Note:  Since  these  subroutines  have  classification  code  LTI-1 
or  is  the  input  and  output  waveforms  are  in  the  fre¬ 
quency  domain. 


4-115 


ivn  SuebOOUM  iHlf-T(X,Y,A,fsJ 

J-yu.  oC'MMOto/bLK  1  /  bn.  1(200 1 

Ivu.<  tOUiVAUNCt  IbKlIltJO)  ,TO>  t (8K1 ( 1 fc U ,TH7 ) ,  IbKl  I  lt2)  .Tjll  I 

lvu-,  *  ,U*1U(S),S1MF0) 

IWj  LilBthiluN  X(  1),  V(  1  >,A«1>  ,B(  1  I 

UATA  NlvbtNlV‘,fNl«»5,Nl*»6/-3.-2.-XtO/ 

iw/  UA1A  Hi  .fcX/o. 2631853,2. VI03fc30t-U/  1780 

Ivuo  o 

ICCN=U 

IWu  I  UK-1  o 

iva  oo  to  to 

lvli  c 


lvli 


tNTKY  fcbHIM  <X»  Yt  A,8  » 


Ob/ll/7!>  1NPU1  L1SI1NI, 

card  nu  ’*** 


AU1UFL0-  CHARI  SCI  -  l-WO/SCL  RAOSIH 
CCMtNTS 


***♦ 


ls>  It 
191t» 
Wlo 
i  w  1  7 
it  it 

It  it 
it*.  U 
It,  1 
it*.  *> 

At,  i 
Lt,w 
It*”? 
i  t  *.  o 
it*/ 

1 1 ,  o 
h," 

it  * i 
it  •»*: 

1 1  s 

it  t 

At 

1‘-  •*' 
it  .  .i 
It  ‘A 
't 

i,t  **o 


U.uN-» 
i.U  Hi  7 

tNlKV  SHlHI><X,Y,A,b> 

UUN  =  l 

1  L  K  -  1  0 

l.U  I>.  lu 

LNlRY  «SII1  lit  X,Y,A,iJl 
U.l  N-  1 
1.1,  U-  *  t 

Hi  .,!•  I  i  -  ll  liuL  »AIU»*»3»  » 

II- (  II  K.  1 1  .O.f  )*>U  1 11  *!u 
NHi.'i 

u  a;.s-  i 

..  -Hi- 1  b.*  *  i 

K  V  t  l  1 

In)  -  I  nl /•>•»*■ 

>>H„  -  ;.  I.-A  ♦  1  HT 

i- 1  l  V  M- -  X(N)  Vj  1*1  L.K 
it  h:.~  <  i  f*  r  f  Hi ) )  *  P I ’ 

i  ,,  y_.:  11.1  LH,-AlNT«lMPHtl*PK’ 

L\  ,  t  L  -  ) 

„M'H-  MNUUL^I 
i  --  v  l  a  i  ► l>  l 
l»i-  jlNU'bl 


i*. ..  i 


i.i.  I> 


\t  o  \  l  1  i'l  iin. 

*  *  *♦ 


1 1  Nit  M  l 


t  Wu/iCL  f  Al.UP 

V*  tfl* 

**** 


1*1  t  -  1  #N*  1  l  £ 

i l nr  =  XI k  ) 

*-«f< i-A<r  J*  J]MMCb-Y(K  >*Sfc 
L(K)-tUhY(U*lUUf!PnN 
1 r hr -l L 

t  i-M.  L  L  “«»’>♦  All.  t  l 

t  l  +lt  Kh  ♦/ Ll  t  t 

K-N*M't  L 
H.  V  tl  i<l  if*C. 

•  «.V.  ]»  l.rK.  -f.1  .Mil  ll  oi  ( 


1  ►  — 

r\  ~ i ,r  I 

:>* 

M*l  l- 

- 1 

r  •  r‘  =.  r  » 

-ULM 

L  an'- 

-‘LllIh 

i  l#  1 1 . 

4  ■' 

i  ■  \i ; 

.  I..;s 

•  l  -  ;s  ( »\  i  **  t 

*  1  .iw 

4-  >-  A(..**  - 

».  <  r  *•* 

-  »-  All.!-. 

(  i  ‘ 

:'l^l  .4*  *• 

l  .1  *■•  i  (  .  i  *  • 
t  .  i  **  .  I  .  /  1  ,  I  *  . 


I 


SUBROUTINE  SWPINT 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

SWPINT  LTI-4  434,435 

NCSWPI  LTI-4  436,437 

2 .  PURPOSE: 

This  subroutine  is  used  to  simulate  a  video  post 
detection  sweep  integrator. 

3.  INPUT  PARAMETERS: 

Name  0/R  T  Description 

FBCK  R  F  Feedback  coefficient 

4.  CALLING  SEQUENCES: 

CALL  SWPINT  (DIN,  DOUT) 

Where:  DIN  contains  the  Input  Waveform 

DOUT  contains  the  Output  Waveform 

The  range  bin  storage  array  (SRI)  is 
cleared  before  execution  begins. 

CALL  NCSWPI  (DIN,  DOUT) 

Where:  DIN  contains  the  Input  Waveform 

DOUT  contains  the  Output  Waveform 

The  range  bin  storage  array  (SRI)  is  not 
cleared  prior  to  execution  except  for  the 
first  execution  of  this  subroutine. 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  maximum  number  of  range  bins  simulated  is  2048. 
If  more  than  2048  range  samples  are  contained  in 
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the  input  array,  only  the  first  2048  are  processed. 


b.  Flow  Chart:  Page  9-204 

c.  Cross  Reference  Table:  Page  9-234 
6.  THEORY  OF  OPERATION 


The  block  diagram  of  the  video  sweep  integrator  simu¬ 
lated  by  this  module  is  shown  in  Figure  SWPINT-1.  The 
Z-plane  transfer  function  for  each  range  bin  is  given  by  the 
following  expression: 

T (Z)  =  Z 

Z  -  FBCK 

The  delay  represented  by  the  Z  operator  is  determined  by 
the  radar  pulse  repetition  interval. 
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SUBROUTINE  TARGET 


1.  MODULE  IDENTIFICATION: 


Name 

Classification  Code 

Reference  Number 

TARGET 

LTI-2 

501 

TGTNCL 

LTI-2S 

502 

2.  PURPOSE: 

This  subroutine  simulates  a  target  represented  by  a 
set  of  discrete  scattering  centers. 


3 .  INPUT  PARAMETERS : 


Name 

F0 

FEXT 

FI 

HTGT 

RTGT0 

ANGTGT 

TORI NT 


O/R  T 

R  F 

R  F 

R  F 

0  F 

0  F 

0  F 

0  F 


Description 

RF  Center  Frequency 

Simulation  bandwidth 

Frequency  increment  of  the 
output  transfer  function 

Height  of  Target 

Range  to  target  (R  in  Figure  1) 
(Used  to  compute  target  return 
variation  as  a  function  of 
range  and  altitude) 

Angle  from  the  radar  reference 
line  to  the  target  (Gin  Figure  1) 

Angle  from  Reference  line  #1  to 
the  target  coordinate  system  ref¬ 
erence  line  (0  in  Figure  1). 
Reference  line  #1  is  parallel 
to  the  radar  coordinate  system 
reference  line  and  passes 
through  the  center  of  target 
rotation. 
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Name  0/R  T  Description 

NSCAT  R  I  Number  of  scatterers 


TSCAT  R  F 


R000  R  F 


TGTVEL  0  F 


Array  containing  the  location 
and  radar  cross  section  of  each 
scatterer.  TSCAT  (l,k)  is  the 
radar  cross  section  of  the  kth 
scatterer(cr,  in  Figure  1). 
TSCAT(2,k)  is  the  range  coordi¬ 
nate  of  the  kth  scatterer  in  the 
target  coordinate  system  (r^  in 
Figure  1).  TSCAT(3,k)  is  the 
angle  coordinate  of  the  kth 
scatterer  in  the  target  coordi¬ 
nate  system  (0-^  in  Figure  1). 

Range  of  target  in  the  simula¬ 
tion  (used  to  compute  time  delay 
only) 

Target  radial  velocity 


TIME  0  F 


Elapsed  time  since  beginning  of 
simulation 


4.  CALLING  SEQUENCES: 

CALL  TARGET  (X,Y) 

Where:  X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

CALL  TGTNCL  (X,Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 


5.  RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  discrete  scattering  center  representation 
used  herein  is  obtained  by  imaging  coherent 
short  pulse  measurements  made  on  the  target  or 
a  model  of  the  target.  The  imaging  procedure 
is  described  in  Appendix  H  of  this  report. 
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h.  The  target  representation  implemented  herein 
includes  the  pulse  to  pulse  phase  shift  due  to 
target  motion  with  respect  to  the  radar  system. 
However,  the  intrapulse  phase  shift  is  not 
included  in  this  model. 

c.  Flow  Chart:  Page  9-109 

d.  Cross  Reference  Table:  Page  9-222 
6.  THEORY  OF  OPERATION 


The  discrete  scatterer  target  scintillation  model 
is  based  on  the  premise  that  a  target  can  be  represented 
as  an  ensemble  of  scattering  centers.  These  scattering 
centers  are  assumed  to  be  sufficiently  small  that  they 
can  be  represented  by  impulses.  Therefore,  this  model  can 
be  described  by  the  following  expression: 


N 


where  o n  is  the  radar  cross  section  of  the  nth 
scatterer 

tn  is  the  time  delay  of  the  nth  scatterer 
<5(  )  is  the  Dirac  delta  function. 

The  Fourier  transform  of  the  above  expression  is  given  by 
the  following  equation: 


n=l 


In  the  target  model  the  following  equations  were  used  in 
evaluating  the  target  response  in  the  frequency  domain. 


4-124 


NSCAT 


X(K)  =  Ge  (ELANG)  *  G^(ANGTGT)*  ^  ^  ° 3  C0S  27rK(  FI+FSTRT)  tn 

J  =  1 
NSCAT 

Y(K)  =  G€  (ELANG)  *  G^(ANGTGT)*  SIN  2jtK(FI+FSTRT)  tn 

J=1 


where:  K  is  the  frequency  index 

FI  is  the  frequency  increment 
FSTRT  is  the  starting  frequency 
FSTRT  =  FO  -  FEXT/2 

G ^  is  the  antenna  azimuth  gain  function 
Gf  is  the  antenna  elevation  gain  function 

ELANG  =  SIN'1  (HTGT/RTGT) 

tn  is  the  time  displacement  of  the 
scatterer 

tn  =  0  119-g^  RTGT  +  TSCAT(2,J)*COS(e-0-TSCAT(3, J)) 
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Figure  TARGET-1  Target  Model  Geometry 
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SUBROUTINE  WEITRE 


1.  MODULE  IDENTIFICATION: 


Name 

WEITRE 

WEITCP 

WEITMP 

2.  PURPOSE: 


Classification  Code 

Reference  Number 

LTI-1  or  LTV-4 

221 

LTI-1  or  LTV-4 

222 

LTI-1  or  LTV-4 

223 

This  subroutine  is  used  to  represent  (1)  a  linear 
time  variant  device  in  terms  of  its  measured  transfer 
functions  or  (2)  a  linear  time  varying  device  in  terms 
of  its  measured  gain  (complex)  as  a  function  of  time. 


3.  INPUT  PARAMETERS: 
Name  O/R  T 

NPWT  R  I 

WT  R  F 


Description 

Number  of  points  used  to  specify 
the  weighting  function 

Array  containing  the  specified 
weighting  function.  The  specifi¬ 
cation  for  the  Jth  sample  is 
the  following: 

WEITRE  or  WEITCP 


WT 

d,J) 

=  Real  component 

WT 

(2,J) 

=  Time 

WT 

(3 ,  J) 

=  Imaginary  component 
(WEITCP  only) 

WEITMP 

WT 

d,J) 

=  Gain 

WT 

(2,J) 

=  Time 

4-129 


WT  (3,J)  =  Phase  angle 


ORIG  0  F  Amount  the  independent  variables 

of  the  weighting  function  is  to 
be  shifted. 

4.  CALLING  SEQUENCES: 

Real  Weighting  Function 
CALL  WEITRE  (X,  Y,  $mmmm) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

Complex  Rectangular  Weighting  Function 
CALL  WEITCP  (X,  Y,  $mmmm) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

Complex  Polar  Weighting  Function 
CALL  WEITMP  (X,  Y,  $mmmm) 

Where:  X  contains  the  Input  Waveform  -  M 

Y  contains  the  Input  Waveform  -  P 

X  contains  the  Output  Waveform  -  M 

Y  contains  the  Output  Waveform  -  P 

$mmmm  is  the  number  of  the  statement  to  which 
control  is  transferred  if  an  error  is  detected, 
in  the  input  data. 

5 .  RESTRICTIONS,  REQUIREMENTS.  AND  MISCELLANEOUS  DATA 

a.  The  input  weighting  function  independent  variable 
must  be  monotonica lly  increasing,  except  that  two 
adjacent  points  may  be  specified  with  the  same 
value.  In  this  case,  if  the  input  value  to  be 
weighted  falls  exactly  on  that  point,  the  first 
dependent  value  specified  will  be  used. 

b.  If  the  extent  of  the  weighting  function  does  not 
include  any  part  of  the  input  data,  the  subroutine 
will  perform  a  nonstandard  exit. 
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c.  If  the  extent  of  the  weighting  function  contains 
only  a  portion  of  the  input  data,  all  other  input 
points  are  set  to  zero. 

d.  The  subroutine  performs  a  linear  interpolation 
between  specified  points. 

e.  Flow  Chart:  Page  9-57 

f.  Cross  Reference  Table:  Page  9-215 
6.  THEORY  OF  OPERATION 


The  basic  mechanization  equations  for  each  entry 
point  are  as  follows: 

WEITCP: 

X(J)  =  X(J)*XWT  -  Y ( J)  *  YWT 
Y  (J)  =  X  (J) *YWT  +  Y  (J)  *  XWT 

A  table  lookup/linear  interpolation  scheme  is  used  to 
determine  the  complex  weight,  XWT  +  j  YWT.  The  value  of 
the  independent  variable  used  in  computing  the  weight  for 
the  Jth  sampie  is  given  by  the  following  expression: 

VIV  =  XORIG  +  J  •  DEL  -  ORIG 

where:  XORIG  =  Independent  variable  for  first  element  of 

the  input  arrays. 

DEL  =  Independent  variable  increment  between 
samples  of  the  input  array. 


WEITRE: 

X(J)  «  X(J)  *  XWT 
Y (J)  -  Y (J)  *  XWT 

The  procedure  for  determination  of  the  weight,  XWT,  is 
the  same  as  that  described  above  except  YWT  is  not  computed. 

WEITMP : 

X(J)  =  X(J)  *  AMPL 
Y  (J)  =  Y(J)  +  PHAS 
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A  table  lookup/interpolation  scheme  is  used  to  determine  the 
complex  weight  AMPL  *  EXP  (PHAS).  The  procedure  for  determin¬ 
ing  value  of  the  independent  variable  used  in  computing  the 
weight  is  described  under  WEITCP. 
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SUBROUTINE  WVGU1D 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

WVGU1D  LTI-1  510 

2.  PURPOSE : 


This  routine  is  used  to  simulate  a  length  of  waveguide 
for  the  purpose  of  predicting  the  effect  of  nonlinear 
phase  behavior  on  waveforms. 


3.  INPUT  PARAMETERS : 
Name  O/R  T 


Description 


RFF  R  F 

CFREQ  R  F 

XWLENG  R  F 

4,  CALLING  SEQUENCES : 

CALL  WVGUID  (X,  Y) 


Center  frequency  of  the  electro¬ 
magnetic  wave  propagating  in  the 
waveguide 

Cutoff  frequency  of  the  waveguide 
Length  of  the  waveguide 


Where: 


X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 


5,  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  References: 

Atwater,  H,  A. ;  Introduction  to  Microwave  Theory. 
McGraw-Hill,  New  York,  162,  Chapter  3, 

b.  Flow  Chart:  Page  9-116 

c„  Cross  Reference  Table:  Page  9-223 
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6.  THEORY  OF  OPERATION 


A  hollow  rectangular  perfectly  conducting  waveguide 
excited  in  only  one  of  its  modes  has  a  propagation  function 
which  is  described  by  the  following  expression: 

/3g(f)  =  ^k2  -  kc2 

where:  k  =  2  n  f  VMoeo 

kc  =  2  tt  f c  V/x0e o 

fc  =  cutoff  frequency 

Substituting  —2^.  =  /i0co  the  following  expression  is  obtained: 

c 

V£>  • 

Therefore  the  transfer  function  for  a  length  of  wave¬ 
guide  considering  only  the  phase  term  is  given  by  the 
following  expressions: 

o(f)  -  e-i  2  V f2  -  fc2  ;  fc  <  f 

-  e  “  Z  Vfc2  -  E2  !  £  S  fc 
where  z  is  the  length  of  the  waveguide. 
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SECTION  5 


CONNECTION  MODULES 


This  section  includes  all  those  modules  that  bind  the 
stimulus/transfer  function  modules  together  in  order  to  form 
the  core  of  the  simulation.  The  functions  performed  by  these 
modules  include  array  addition  or  array  multiplication  (CONV) 
and  conversion  between  time  and  frequency  domain  representa* 
tions  (ZFFT) . 


SUBROUTINE  CONV 


1.  MODULE  IDENTIFICATION: 


Name 

Classification  Code 

Reference  Number 

CONV 

Connection 

204 

CONVMP 

Connection 

205 

DIVA 

Connection 

206 

ADDA 

Connection 

207,239 

2.  PURPOSE: 

This  subroutine  serves  to  connect  simulation  modules 
together  by  performing  array  addition,  multiplication  and 
division . 


3.  LNPUT  PARAMETERS: 
None 


4.  CALLING  SEQUENCES: 

Complex  Rectangular  Array  Multiplication 


CALL  CONV  (XI,  XQ,  FI,  FQ ,  $mmmm) 


Where : 


XI  contains  the  Input  Waveform  -  R 
XQ  contains  the  Input  Waveform  -  I 
FI  contains  the  Input  Waveform  -  R 
FQ  contains  the  Input  Waveform  -  I 
XI  contains  the  Output  Waveform  -  R 
XQ  contains  the  Output  Waveform  -  I 


Complex 
Array  #1 
Complex 
Array  #2 


mmmm  is  the  return  address  for  a  non-standard 
re  turn 


XI(J)  =  XI(J)  *  FI( J)  -  XQ( J)  *  FQ( J) 
XQ(  J)  =  XI(J)  *  FQ( J)  +  XQ( J)  *  FI( J) 


Complex  Polar  Array  Multiplication 
CALL  CONVMP  (XI,  XQ,  FI,  FQ,  $mmmm) 

Where:  XI  contains  the  Input  Waveform  -  M  Complex 

XQ  contains  the  Input  Waveform  -  P  Array  #1 

FI  contains  the  Input  Waveform  -  M  Complex 

FQ  contains  the  Input  Waveform  -  P  Array  #2 

XI  contains  the  Output  Waveform  -  M 

XQ  contains  the  Output  Waveform  -  P 

mmmm  is  the  return  address  for  a  non¬ 
standard  return. 

XI(J)  =  XI(J)  *  FI ( J) 

XQ(J)  =  XQ( J )  +  FQ( J) 

Complex  Polar  Array  Division 

CALL  DIVA  (XI,  XQ,  FI,  FQ ,  $mmmm) 

Where:  XI  contains  the  Input  Waveform  -  M  Complex 

XQ  contains  the  Input  Waveform  -  P  Array  //I 

FI  contains  the  Input  Waveform  -  M  Complex 

FQ  contains  the  Input  Waveform  -  P  Array  #2 

XI  contains  the  Output  Waveform  -  M 

XQ  contains  the  Output  Waveform  -  P 

mmmm  is  the  return  address  for  a  non-standard 
re  turn 

XI(J)  =  XI ( J) /FI ( J) 

XQ(J)  =  XQ( J) -FQ( J) 

Array  Addition 

CALL  ADDA  (XI,  FI,  $mmmm) 

Where:  XI  contains  the  Input  Waveform 

FI  contains  the  Input  Waveform 

XI  contains  the  Output  Waveform 

mmmm  is  the  return  address  for  a  non-standard 
return 

XI (J)  =  XI (J)  +  FI( J) 
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5. 


RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 


a.  The  difference  in  the  independent  variable  incre¬ 
ment  of  the  input  arrays  must  be  less  than  or 
equal  to  10“^. 

b.  The  output  is  set  to  zero  for  any  cases  where  the 
independent  variables  of  the  two  arrays  do  not 
overlap . 

c.  When  polar  operations  are  performed  the  phase  angle 
is  measured  in  degrees. 

d.  If  the  arrays  do  not  overlap,  or  if  they  must  be 
shifted  more  than  4097  increments,  the  problem 
will  abort  and  the  non-standard  return  will  be 
used . 

e.  Flow  Chart:  Page  9-60 

f.  Cross  Reference  Table:  Page  9-216. 

6.  THEORY  OF  OPERATION 


Depending  upon  the  entry  point  through  which  the  sub¬ 
routine  is  called,  a  control  flag  is  set  to  an  integer 
number  1  to  4.  This  number  indicates  to  the  subroutine  the 
type  of  operation  which  is  to  be  performed,  either  rectan¬ 
gular  multiplication,  polar  multiplication,  polar  division, 
or  addition.  After  the  flag  is  set,  the  two  arrays  are 
examined  (beginning  at  statement  number  1000)  to  verify  that 
the  values  of  the  increments  of  the  tw •>  arrays  are  within 
1/10000  of  each  other.  The  arrays  are  inspected  to  verify 
that  they  overlap,  at  least  in  part,  and  one  array  is  shifted, 
if  required,  to  align  the  two  origins.  If  the  .  ncrements 
are  not  within  1/10000  of  each  other,  or  if  the  arrays  do 
not  over  Lap,  or  if  the  origins  of  the  two  arrays  are  not 
within  4097  increments  of  each  other,  the  problem  is  aborted. 
Otherwise,  the  proper  arithmetic  manipulation  is  performed, 
according  to  the  value  of  the  control  flag  and  the  results 
are  placed  in  the  X  array(s).  The  F  array(s)  remain  unmod¬ 
ified  after  exit  from  the  subroutine.  For  the  divide  opera¬ 
tion,  the  X  complex  array  is  the  dividend  and  the  F  complex 
array  is  the  divisor. 
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SUBROUTINE  ZFFT 


1.  MODULE  IDENTIFICATION: 


Name 


Classification  Code  Reference  Number 


ZFFT 

ZIFFT 

IFT 


Connection 

Connection 

Connection 


202 

203 

Not  User  Referenced 


2.  PURPOSE: 

This  subroutine  performs  the  finite  discrete  Fourier 
Transform  (ZFFT)  and  inverse  Fourier  Transform  (ZIFFT  or 
IFT)  of  a  sequence  of  input  data  samples.  The  mechaniza¬ 
tion  is  based  on  the  Fast  Fourier  Transform  (FFT)  algo¬ 
rithm  developed  by  Langdon  and  Sande  from  the  approach 
of  J.  W.  Tukey  and  J.  Cooley. 


3.  INPUT  PARAMETERS: 

Name  0/R  T  Description 

N2  R  I  Power  of  2  which  determines  the  total 

number  of  prints  (NTHPOW)  transformed 
by  the  FFT. 


FSHIFT  R  F  Frequency  shift  to  be  applied  to  fre¬ 
quency  domain  representation  of  the 
waveform. 


ICPLXI  0  I  Control  flag  which  indicates  the 

nature  of  the  inverse  transform  (ZIFFT) 
input  data 

=  1  complex  waveform  Z  =  X  +  jY 
=  0  real  waveform  Z  =*  X 


ICPLXF  0  I  Control  flag  which  indicates  the  nature 

of  the  forward  transform  (ZFFT)  input 
data 

=  1  complex  waveform  Z  *  X  +  jY 
■  0  real  waveform  Z  *  X 

SIMBW  R  F  Width  of  output  spectrum  when  forward 

transform  (ZFFT)  is  performed 
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Name  0/R  T  Description 

INORM  R  I  Normalization  flag 


2FFT  ZIFFT 


»  0; 

"ti 

r 

1  N 

=  FI 

=  1; 

fN  " 

1.0/NP 

fN 

*  FI 

=  2; 

fN  - 

TI 

fN 

1  *0/ (2**N2) 

=  3; 

fN  - 

1 *°/NP 

fN 

=  1 *0/ (2**N2) 

-  4; 

fN  = 

1.0 

fN 

-  1.0 

4.  CALLING  SEQUENCES: 

Fourier  Transform 
CALL  ZFFT  (X,Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 
Increase  Fourier  Transform 

CALL  ZIFFT  (X,  Y) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

Inverse  FFT  only 
CALL  I FT  (X,  Y) 


l 
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Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 

X  contains  the  Output  Waveform  -  R 

Y  contains  the  Output  Waveform  -  I 

This  entry  point  is  used  in  computing  antenna  patterns 

and  is  called  by  ANTPAT  module  only. 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  maximum  value  of  N2  is  13  which  gives  8192 
transformed  samples. 

b.  The  number  of  points  transformed  by  the  FFT  algo¬ 
rithm  is  2  raised  to  an  integer  power.  If  the  num¬ 
ber  of  input  samples  is  not  equal  to  a  power  of  2, 
then  zeros  are  used  to  fill  in  the  array  elements 
not  defined  by  the  input  data;  i.e.,  zeros  are 
placed  in  the  input  array  locations  from  NP  +  1 

to  2**N2 . 

c.  When  the  Forward  Transform  is  executed  the  output 
spectrum  is  shifted  by  FSHFT  and  has  a  width  of 
SIMBW:  i.e.,  the  output  spectrum  extends  from  FSHIFT  - 
SIMBW/2  to  FSHIFT  +  SIMBW/2. 

d.  When  the  Forward  Transform  is  executed  the  origin 
of  the  independent  variable  X  (N194)  is  set  to  the 
discrete  frequency  nearest  FSHIFT  -  SIMBW/2  since 
the  FFT  output  consists  of  a  discrete  set  of  fre¬ 
quency  samples.  When  the  Inverse  Transform  is  exe¬ 
cuted  the  origin  of  the  independent  variable  X(N194) 
is  set  to  zero. 

e.  The  independent  variable  increment  is  set  to  1.0/ 
(X(N195)*(2**N2POW) )  for  both  the  forward  and  inverse 
transform. 

f.  When  the  inverse  transform  is  executed  the  input 
spectrum  is  shifted  by  FSHIFT  before  the  FFT  is 
executed.  The  value  of  FSHIFT  when  the  inverse 
transform  is  performed  can  be  different  from  the 
value  used  when  the  forward  transform  was  performed. 


5-9 


The  following  relationship  must  be  satisfied 
-  FS/2<FSHIFT  +  SIMBW/2<  FS/2 

g.  Flow  Chart:  Page  9-125 

h.  Cross  Reference  Table:  9-224 
THEORY  OF  OPERATION 


a.  ZFFT  -  This  entry  point  connects  the  time  domain 

representation  of  a  waveform  into  the  corresponding 
frequency  domain  representation.  The  Fourier  trans¬ 
form  relationship  is  given  by  the  following  expression 


S(f) 


s(t) 


-1 2r f t 
e 


dt 


The  basic  mechanization  equation  used  in  computing 
the  discrete  Fourier  transform  is  the  following: 


NP-1 


0 


where:  k  is  the  frequency  index  0  k  NTHPOW- 
i  is  the  time  index  0  i  NP-1 
f  is  the  normalization  factor 
W  «  exp  ( j  2  VNTHPOW) 

In  general  s^  and  are  complex  arrays. 

b.  ZIFFT  and  IFT  -  These  entry  points  connect  the  fre¬ 
quency  domain  representation  of  a  waveform  into  the 
corresponding  time  domain  representation.  The 
inverse  Fourier  transform  relationship  is  given  by 
the  following  expression: 


s(t) 


S(f)  e 


j  2*f  t 


dt 
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The  basic  mechanization  equation  used  in  computing 
the  inverse  discrete  Fourier  transform  is  the 
following: 

NTHPOW-1 

ki 

0  <  i  <  NTHPOW-1 

0 

The  only  difference  between  ZIFFT  and  IFT  is  that 
IFT  performs  the  inverse  FFT  only  and  no  normaliza¬ 
tion  is  performed.  The  FFT  mechanization  utilized 
in  this  module  evaluates  the  above  expression  directly  . 
In  order  to  compute  the  forward  transform  using  the 
same  inverse  FFT  the  following  expression  is  used. 


s.  w 

k 
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SECTION  6 


PERIPHERAL  MODULES 


This  section  includes  all  those  modules  that  do  not 
participate  in  the  actual  simulation  process,  but  serve  as 
the  interface  between  the  simulation  and  the  user.  In  some 
cases  this  is  as  simple  as  outputting  data  to  a  disc  file. 

In  other  cases  the  data  is  processed  to  extract  some  para¬ 
meter,  e.g.  energy  in  a  waveform,  or  to  map  data  into  a 
different  representation  space,  e.g.  the  generation  of  a 
probability  histogram.  The  following  modules  are  included 
in  this  group: 

CLINT 

CUMDIS 

ERGYCP 

PLOTTR 

PLTFMT 

PTLIST 

RTOPDB 

SCANNR 

Also  included  in  this  group  is  the  module  SPCAVG,  which 
is  located  in  Volume  III. 

The  following  peripheral  modules  are  related  to  the 
bi static  target  model  and  target  imaging  and  are  included 
in  Volume  I,  Part  3: 

ERRPRC 

(BISTATIC  ANTENNA  INITIALIZATION) 
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SUBROUTINE  CLINT 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

CLINT  PERIP.  302 

2.  PURPOSE: 

This  subroutine  initializes  the  clutter  model 
scatterers  and  parameters. 

3.  INPUT  PARAMETERS: 

Name  O/R  T  Description 

LAMBDA  0  F  Wavelength  of  the  transmitted 

waveform.  This  parameter  is 
used  to  compute  the  doppler 
shift  of  weather  or  chaff 
c lutter . 

TCELL  R  F  Range  dimension  of  clutter  element. 

IDMP  0  I  Output  data  print  flag. 

IDMP  *  0  (a)  Scatterer  arrays  are 
printed  (CLINT) 

(b)  Antenna  Azimuth  and 
Elevation  arrays  are 
printed  (cluttr) 

IDMP  =  1  Clutter  impluse 

response  is  printed 
by  (cluttr). 


NRCS 

R 

I 

Scatterer  RCS  probability  dis¬ 
tribution  flag. 

RWPH 

0 

F 

Scatterer  random  walk  maximum 
phase  shift.  Represents  the 
random  motion  of  cluttr  from  one 
transmit  cycle  to  the  next. 

WNDVEL 

0 

F 

Speed  of  the  clutter  volume  due 
to  wind,  meters/second. 
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p 


Name 

0/R 

T 

Description 

TI 

R 

F 

Time  increment  between  samples 
(used  to  compute  number  of 
storage  cells  required  by  clutter 
impulse  response) . 

VELANG 

0 

F 

Relative  angle  of  clutter  velocity 
vector  with  respect  to  the  radar 
site. 

RNEXT 

0 

F 

Range  dimension  of  the  clutter 
volume  in  monoseconds. 

RN000 

0 

F 

Clutter  volume  starting  range. 

AZEXT 

0 

F 

Azimuth  dimension  of  the  clutter 
volume  in  degrees 

AZ000 

0 

F 

Clutter  volume  starting  azimuth 
angle . 

MM 

0 

I 

Number  of  clutter  volume  azimuth 
increments. 

ELEXT 

0 

F 

Elevation  dimension  of  the  clut¬ 
ter  volume  in  degrees. 

EL000 

0 

F 

Clutter  volume  starting  elevation 
angle. 

NN 

0 

I 

Number  of  clutter  volume  eleva¬ 
tion  increments. 

4.  CALLING  SEQUENCES: 

CALL  CLINT  ($NNNN) 

WHERE:  NNNN  Statement  number  for  abnormal 

return 


*  . 


( 
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* 


5.  RESTRICTION.  REQUIREMENTS.  MISCELLANEOUS  DATA: 

a.  The  clutter  model  was  split  into  two  subroutines 
in  order  that  the  initialization  process  would 
be  separate  and  easily  changed  to  allow  modeling 
of  more  complex  clutter  environments. 

b.  The  total  number  of  scatterers  (NSCAT)  is  limited 
only  by  the  amount  of  available  disc  space  and  is 
approximately  equal  to  160  times  the  number  of 
blocks  of  available  disc  space.  The  number  of 
scatterers  is  equal  to  NN  times  MM  times  KK,  where 
KK  -  RNEXT/TCELL.  Other  than  their  product;  NN, 

MM  and  KK  have  no  limitations. 

c.  The  random  number  generator  function  (RRAND) 
is  called  from  this  subroutine. 

d.  The  clutter  scatterer  element  arrays  are  stored 
in  a  random  disc  file,  FC  *  02,  the  record  size 
is  500  words. 

e.  Flow  Chart:  Page  9-198 

f.  Cross  Reference  Table:  Page  9-234 

6.  THEORY  OF  OPERATION 


This  program  can  be  divided  into  two  functional  sec¬ 
tions  and  each  will  be  discussed  in  the  following  paragraph. 

a.  User  parameter  testing  and  clutter  model  parameter 
calculation.  Refer  to  Figures  CLINT  -  1(a)  and  1(b). 


XVLANG  The  angle  between  the  wind  velocity  vector 
and  a  radial  line  from  the  first  clutter 
scatterer  element  (0,0,0)  to  the  radar. 

XVLANG  =  (AZ000  -  VELANG) /57 .29578 

DOPFRQ  The  maximum  doppler  frequency  shift  due  to 
the  wind. 

DOPGRQ  =  -2.0  *  WNDV EL/ LAMBDA 

If  LAMBDA  is  zero  this  calculation  is  not 

performed  and  DOPFRQ  is  set  to  zero. 
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RADAR  SITE 


Figure  CLINT- lb  CLUTTER  MODEL  GEOMETRY 

O'lANE  VIEW) 


h-b 


MM 

NN 

KCELL 

KK 

NBLKS 

DELAZ 

DELEL 

ICFLG 

MODE 


If  this  parameter  was  not  specified  it  is 
set  to  1. 

If  this  parameter  was  not  specified  it  is 
set  to  1. 

The  total  number  of  storage  cells  required 
to  represent  the  clutter  impulse  response. 

If  this  parameter  exceeds  4096  (half  the 
available  storage)  then  RNEXT  is  auto¬ 
matically  reduced  to  make  KCELL  <4096. 

KCELL  =  RNEXT /TI 

Number  of  clutter  volume  range  increments. 

KK  =  RNEXT/TCELL 

The  number  of  random  disc  records  required 
to  store  the  clutter  scatterer.  This 
parameter  is  arbitrarily  limited  to  200  and 
can  be  increased  providing  the  control 
cards  for  DSRN  02  are  changed  accordingly. 

Azimuth  increment  between  clutter  scatterers 

DELAZ  =  AZEXT/MM 

Elevation  increment  between  clutter 
scatterers . 

DELEL  =  ELEXT/NN 

This  parameter  is  set  to  1  if  the  clutter 
model  has  been  successfully  initialized. 

This  is  tested  before  CLUTTR  can  be  executed 

This  parameter  is  used  to  indicate  if 
stationary  (MODES)  or  time  varying  (M0DE=2) 
clutter  is  to  be  modeled.  When  stationary 
clutter  is  modeled  certain  simplifications 
can  be  made  in  the  processing  performed  in 
the  cultter  model  routine  CLUTTR. 
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b.  Each  clutter  scatter  element  Is  represented 

by  a  radar  cross  section  value  and  a  phase  angle. 
These  two  parameters  are  converted  to  rectanqular 
representation  and  stored  in  the  arrays  CLUX  and 
CLUY. 

CLUX(J)  -  A*RRAND(8) 

CLUY ( J)  -  A*DUM 

where  A  -  v/rCS(J)  »  /RRAND(NRCS) 

RRAND(8)  «  COS (0) 

DUM  *  SIN(0) 

6  is  a  sample  of  a  random  uniform 
distribution  3OO.O<0£O.O 
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SUBROUTINE  CUMDIS 


1.  MODULE  IDENTIFICATION: 


Name 


Classification  Code  Reference  Number 


CUMDIS 

CUM2 

OUTCUM 

PDF 

2 .  PURPOSE : 


Peripheral 

Peripheral 

Peripheral 

Peripheral 


208,209 

Not  User  Referenced 
210,211 
212,213 


This  subroutine  processes  a  sequence  of  input  samples 
in  order  to  generate  the  probability  density  histogram 
and  cumulative  probability  distribution. 


3. 


INPUT  PARAMETERS 


Name 

0/R 

T 

TLIM 

R 

F 

BLIM 

R 

F 

NIXF 

R 

I 

INPUT 

PARAMETERS 

Name 

0/R 

1 

NCPACK 

R 

I 

INPUT 

PARAMETERS 

Name 

0/R 

1 

NDPACK 

R 

I 

(CUMDIS): 


Description 

The  upper  limit  of  the  histogram 
independent  variable 

The  lower  limit  of  the  histogram 
independent  variable 

The  number  of  elements  in  the 
histogram 

(OUTCUM): 


Description 

The  number  of  histogram  elements  com¬ 
bined  to  form  one  element  of  the  out¬ 
put  cumulative  distribution 

(PDF): 


Description 

The  number  of  histogram  elements  com¬ 
bined  to  form  one  element  of  the  out¬ 
put  probability  density. 
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4.  CALLING  SEQUENCES: 


Initial  histogram  computation 
CALL  CUMDIS  (DATAIN,  XF) 

Where:  DATAIN  contains  the  Input  Waveform 

XF  contains  the  Output  Histogram 

Add  additional  data  to  histogram 
CALL  CUM2  (DATAIN,  XF) 

Where:  DATAIN  contains  the  Input  Waveform 

XF  contains  the  Output  Histogram 

Cumulative  distribution 
CALL  OUTCUM  (DATOUT,  XF) 

Where:  XF  contains  the  Input  Histogram 

DATOUT  contains  the  cumulative 
distribution 

Probability  density 
CALL  PDF  (DATOUT,  XF) 

Where:  XF  contains  the  Input  Histogram 

DATOUT  contains  the  probability  density 
5.  RESTRICTIONS .  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  The  number  of  points  NIXF  in  the  histogram  XF 
cannot  exceed  8000. 

b.  Any  value  of  the  input  array  which  lies  below  the 
bound  of  BLIM  will  be  placed  in  the  lowest  element 
of  the  histogram.  Any  value  of  the  input  array 
which  lies  above  the  bound  of  TLIM  will  be  placed 
in  the  highest  element  of  the  histogram. 
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c.  The  variables  NCPACK  and  NDPACK  must  be  greater 
than  or  equal  to  1. 

d.  Flow  Chart:  Page  9-152 

e.  Cross  Reference  Table:  Page  9-228 
6.  THEORY  OF  OPERATION 

a.  The  width  of  the  histogram  elements  is  determined 
using  the  following  expression 

DINC  =  (TLIM  -  BLIM)  /  NIXF 

b.  The  value  of  each  sample  in  the  input  array  DATAIN 
is  divided  by  DINC  to  determine  the  corresponding 
histogram  element.  The  number  of  samples  which  lie 
within  the  bounds  of  that  histogram  element  is  incre¬ 
mented  by  one. 

c.  After  the  histogram  is  generated  an  entry  through 
OUTCUM  causes  the  cumulative  distribution  to  be 
calculated  and  placed  in  array  DATOUT.  The  number 
of  histogram  elements  combined  to  generate  one  out¬ 
put  sample  is  determined  by  the  variable  NCPACK. 

d.  If  the  density  function  is  desired,  entry  is  made 
through  the  entry  point  PDF.  The  number  of  histo¬ 
gram  elements  combined  to  generate  one  output 
sample  is  determined  by  the  variable  NDPACK. 
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SUBROUTINE  ERGYCP 


MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


ERGYCP 

ERCYRE 

2.  PURPOSE: 


Peripheral 
Periphera 1 


118 

116,117 


This  subroutine  computes  the  energy  contained  in 
a  waveform. 

INPUT  PARAMETERS: 

None 

CALLING  SEQUENCES: 

Complex  waveform 
CALL  ERGYCP  (X,Y) 

Where :  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
NPTS 

Energy  =  X(J)  ‘"''(J)  +  Y(J)  ,f  Y(J) 

J  =  ] 

NPTS  =  Number  of  waveform  samples 

A  =  Independent  variable  increment  between 
samp les 

Real  Waveform 
CALL  ERGYRE  (X) 


Where : 


Energv  = 


X  contains  the  Input  Waveform 
NPTS 

A  X(J)  *  X(J) 

J  =  1 
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a.  The  computed  energy  is  printed  on  the  output  listing. 

The  measurement  unit  is  watt-nonoseconds ,  normally. 

b.  Flow  Chart:  Page  9-70 

c.  Cross  Reference  Table:  Page  9-217. 
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SUBROUTINE  PLOTTR 


I.  MODULE  IDENTIFICATION: 


Name  Classification  Code  Reference  Number 

PLOTTR  PERIP  307,308,309,310 

2.  PURPOSE: 


This  subroutine  converts  the  data  contained  in  the 
input  array  into  a  form  suitable  for  plotting  by  a  HP9820 
desk  calculator.  In  addition,  the  output  can  also  be  listed 
on  a  TTY  via  TSS  or  punched  on  cards. 


3.  INPUT  PARAMETERS: 


Name  0/R  T 

SI  VST  R  F 

XIVRNG  R  F 

LOGFLG  R  I 


Description 

Value  of  the  independent  variable 
at  the  start  of  the  plot 

Range  of  the  independent  variable 
which  is  to  be  plotted 

Control  index  which  indicates  the 
nature  of  the  input  data 


NSKP  0  I  This  variable  controls  the  number 

of  input  samples  shipped  between 
each  plotted  point.  For  example, 
if  NSKP=0  or  1,  each  point  is 
plotted;  if  NSKP=2,  every  other 
point  is  plotted. 

NAUTO  0  I  Control  flag  for  automatic  data 

scanner  which  determines  the  maxi¬ 
mum  and  minimum  values  of  the  de¬ 
pendent  variable.  If  NAUT0=0,  the 
scanner  determines  the  value  of  TH 
and  TL.  If  NAUT0=1,  the  use  pro¬ 
vided  input  parameters  TH  and  TL 
are  used. 


I'H  C  F  Maximum  dependent  variable  value 

to  be  plotted 


6-18 


Name 


Description 


0/R  T 
IC  R  I 

VIL  0  C 

VDL  0  C 

GLBL  0  C 

CALLING  SEQUENCES: 
CALL  PLOTTR  (DV) 
Where:  DV  conta 


Logical  unit  designator  for  output 
file 

Independent  variable  label 

Dependent  variable  label 

Plot  title  and  miscellaneous  data 


ins  the  Input  Waveform 


5.  RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  If  an  output  file  code  (IC)  is  not  specified 
then  no  data  will  be  transferred  and  control 
will  immediately  be  transferred  back  to  the 
calling  program. 

b.  If  the  user  specified  parameters,  XIVST  and 
SIVRNG,  specify  a  plot  that  is  outside  the 
range  of  the  input  data  array  then  the  complete 
input  array  is  processed  and  transferred  to  the 
output  file. 

c.  If  the  file  code  is  negative,  specifying  HP 
paper  tape  format,  then  all  plot  labels  and 
miscellaneous  data  are  suppressed.  Only  TH,  TL. 
XIVOR,  DELTA,  and  NPTS  and  the  plot  data  are 
transferred . 

d.  All  floating  point  output  values  are  in  E13.5 
forma  t . 

e.  Flow  Chart:  Page  9-90 

f.  Cross  Reference  Table:  Page  9-219. 
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6. 


THEORY  OF  OPERATION 


i.  I  pon  entering  the  subroutine  l  ho  paramo  tor  10  is 
tested  to  determine  that  1 0  is  not  zero  and  what 
operating  mode  is  to  bo  used.  If  I C  0  then  the 
HP  paper  tape  mode  (ICON  *  2)  is  selected.  If  IC  0 
then  the  card  image  mode  (IC.0N=1)  is  selected. 

b.  if  (IC0N=1)  the  following  lines  (records)  are 
transmitted  to  the  output  file. 

1 .  Header  line  "**  Punch  Card  Data  Output" 

2.  XT  VST ,  XIVRNG,  XIVFR  (not  used),  LOGFLG , 

VI L  (first  12  characters) 

3.  V'L  (last  60  characters) 

4.  VDL  (first  60  chara  ters) 

5.  VDL  (last  12  characters),  GLBL  (first  48 

characters ) 

c.  The  input  parameters,  XVST  and  XIVRNG,  are  tested 
to  verify  that  the  range  of  the  desired  plot  is 
within  the  limit  of  the  input  array,  DV.  If  the 
desired  plot  extends  past  either  end  of  DV,  then 

t  lie  values  of  XTVST  and  XIVRNG  are  changed  to 
he  compatible  with  the  available  data.  If  the 
desired  plot  is  totally  outside  the  available 
data,  then  the  complete  input  array  is  processed. 

d.  If  NAl’TO  =  0  the  input  data  is  scanned  to  deter¬ 
mine  the  maximum,  TH,  and  the  minimum,  TL,  values 
of  the  dependent  variable  to  be  plotted. 

e.  1F(IC0N=1)  the  following  lines  (records)  are 
transmitted  to  the  output  file. 

1.  GLBL  (last  10  characters) 

2.  NSKP ,  NAl'TO ,  IC 

f.  The  following  lines  (records)  are  transmitted  to 
the  output  file. 

1.  TH,  TL,  XIVOR,  DELTA,  NPTS 

2.  5  dependent  variable  values 
!> .  5  dependent  variable  values 


X-  Last  line  of  plot  data 
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SUBROUTINE  PLTFMT 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

PLTFMT  PERIP  113 

2.  PURPOSE: 

This  subroutine  is  used  to  transform  input  data  which 
is  provided  at  equal  increments  of  sine  0  into  an  output 
array  which  is  equally  incremented  in  9.  The  new  points 
are  determined  using  a  parabolic  interpolation.  In  addi¬ 
tion,  the  input  array  is  tested  to  determine  if  all  input 
points  are  in  visible  space  (independent  variable  between 
+  1),  and  superflous  data  is  eliminated. 

3.  INPUT  PARAMETERS: 


Name 

911 

T 

Description 

DELTHE 

R 

F 

Output  array  independent  variable 
increment 

STARPT 

0 

F 

Lower  limit  of  the  output  array 

STOPT 

0 

F 

Upper  limit  of  the  output  array 

4.  CALLING  SEQUENCES: 

CALL  PLTFMT  (GXIN,  GYIN,  GOUT,  $nnnn) 

Where:  GXIN  contains  the  Input  E-field  -  Real 

Component 

GYIN  contains  the  Input  E- field  -  Imaginary 

Component 

GOUT  contains  the  Output  gain 

nnnn  is  the  statement  number  to  which 

control  is  transferred  when  a  non- 
recoverable  error  is  detected. 

5.  RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  lower  and  upper  limits  of  the  output  array 
must  be  specified  in  ascending  order  in  order  to 
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be  e  f  fee  t  i  ve  .  If  the-  upper  limit  is  less  than  or 
equal  to  the  lower  limit  or  no  limits  are  speci¬ 
fied,  all  valid  input  data  will  be  interpolated. 

b.  A  double  parabolic  interpolation  is  performed  to 
determine  each  output  sample  except  in  the  vicinity 
of  the  end  points. 

c.  Flow  Chart:  Page  9-175 

d.  Cross  Reference  Table:  Page  9-231. 

6.  THEORY  OF  OPERATION 

PLTFMT  is  intended  primarily  for  converting  antenna 
pattern  data  from  equal  increments  of  sine  0,  generated 
by  the  Fast  Fourier  Transform,  to  equal  increments  of  0. 

This  conversion  is  important  in  that  it  allows  the  pattern 
array  to  be  stored  as  an  array  of  gain  points,  equally 
spaced  in  angle.  This  eliminates  the  necessity  of  storing 
a  gain-angle  data  pair,  since  the  location  of  any  point  in 
angle  can  be  determined  from  the  ordinal  number  of  the 
point,  the  origin  of  the  array,  and  the  angle  increment. 

It  is  especially  important  in  converting  antenna  pattern 
data  for  plots,  since  it  is  desirable  that  the  points  be 
equispaced  in  the  independent  variable  direction.  An  addi¬ 
tional  feature  of  PLTFMT  is  that  it  will  automatically 
eliminate  the  extraneous  output  from  the  Fast  Fourier  Trans¬ 
form,  which  is  contained  in  the  region  commonly  known  as 
"invisible  space",  i.e.,  that  region  where  the  dependent 
variable  lies  outside  -1  5  sine  05+1.  Furthermore,  by 
specifying  angle  limits  for  the  region  of  interest,  super- 
flous  data  is  eliminated,  thereby  reducing  processing  time 
and  storage  requirements. 

a.  In  the  first  section  of  the  program  the  input  array 
input  parameters  are  initialized.  The  end  points 
specified  for  the  input  data  are  examined  to  deter¬ 
mine  if  they  are  realistic  and  within  the  boundaries 
of  the  input  array.  If  not,  one  of  three  things 
happens .  (1)  If  the  specified  end  points  exceed 

the  array  limits  or  +  90  degrees,  they  are  reduced 
to  a  useable  value.  (2)  If  the  limits  are  equal 
or  reversed,  i.e.,  start  point  equals  stop  point 
or  start  point  is  greater  than  stop  point,  the 
limits  are  set  to  +  90  degrees  or  to  be  array 
limits,  whichever  is  less.  (3)  If  the  specified 
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end  points  are  in  the  proper  order,  but  do  not 
encompass  any  of  the  array  data,  the  problem  is 
terminated.  If  the  adjusted  end  points  do  not 
encompass  all  of  the  input  array,  the  portion  of 
the  input  array  to  simplify  processing  and  the 
array  definition  parameters  are  adjusted  accordingly. 

b.  The  output  array  elements  are  computed  by  inter¬ 
polating  between  the  input  array  points.  A  double 
parabolic  interpolation  is  used  except  in  the  end 
segments  where  a  single  parabolic  interpolation 
is  used. 

The  double  parabolic  interpolation  is  performed  by 
generating  two  quadratic  equations,  using  three 
adjacent  points  to  generate  each  equation.  Assum¬ 
ing  the  points  are  labeled  XI,  X2,  X3 ,  and  X4,  and 
the  required  location  is  between  X2  and  X3,  then  the 
equations  are  generated  using  XI,  X2,  X3,  and  X2, 

X3 ,  X4.  Once  the  quadratic  coefficients  have  been 
determined,  the  value  of  each  quadratic  at  the 
point  of  interest  is  determined,  and  the  two  values 
are  averaged  to  find  the  value  of  the  dependent 
variable . 

Within  the  end  segments,  adjacent  points  are  not 
available  to  obtain  two  quadratic  equations.  There¬ 
fore,  a  single  quadratic  is  generated  and  the 
necessary  points  are  interp  lated  from  this  single 
equation.  In  all  cases,  if  the  point  of  interest 
falls  exactly  on  a  point  in  the  input  array,  the 
actual  value  of  that  point  is  used. 
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SUBROUTINE  PTLIST 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

PTLIST  PERIP  303,304,305,306 

2 .  PURPOSE : 


This  subroutine  converts  the  data  contained  in  the  input 
array  into  a  compact  form  for  transmission  to  a  remote 
processing  station  which  has  a  CRT  plotter. 


3.  INPUT  PARAMETERS: 
Name  O/R  T 

ST  R  F 

RNG  R  F 

LF  R  I 

NSKP  0  I 


NAUTO  0  I 


TH  OF 


Description 

Value  of  the  independent  variable 
at  the  start  of  the  plot 

Range  of  the  independent  variable 
which  is  to  be  plotted 

Control  index  which  indicates  the 
nature  of  the  input  data 

This  variable  controls  the  number 
of  input  samples  skipped  between 
each  plotted  point.  For  example, 
if  NSKP  *  0  or  1,  each  point  is 
plotted;  if  NSKP  *  2,  every  other 
point  is  plotted. 

Control  flag  for  automatic  data 
scanner  which  determines  the 
maximum  and  minimum  values  of  the 
dependent  variable.  If  NAUTO  ■  0 
the  scanner  determines  the  value 
of  TH  and  TL.  If  NAUTO  -  1,  the 
user  provided  input  parameters 
TH  and  TL  are  used 

Maximum  dependent  variable  value 
to  be  plotted 
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AD-A102  278 


GENERAL  DYNAMICS  FORT  WORTH  TEX  CONVAIR  AEROSPACE  DIV  F/G  17/9 

ENDO  ATMOSPHERIC-EXO  ATMOSPHERIC  RADAR  MODELING.  VOLUME  II.  PAR— ETC (U> 
JUN  76  R  J  HANCOCK  t  F  H  CLEVELAND  F30602-73-C-0380 

RA0C-TR-76-186-V0L-2-PT-1  NL 


UNCLASSIFIED 


Name 

0/R  T 

Description 

TL 

0 

P 

Minimum  dependent  variable  value 
to  be  plotted 

IFCODE 

R 

I 

Logical  unit  designator  for 
output  file 

h.  CALLING  SEQUENCES: 

CALL  PTLIST  (DV) 

Where:  DV  contains  the  Input  Waveform 


RESTRICTIONS.  REQUIREMENTS .  MISCELLANEOUS  DATA 

a.  A  number  of  tests  are  performed  on  the  input  data 

to  ensure  that  a  plot  can  be  successfully  completed. 
If  all  the  necessary  conditions  are  not  satisfied, 
a  condition  code  index  is  set  and  the  plotter  is 
bypassed  with  control  returned  to  the  calling  pro¬ 
gram.  The  following  is  a  list  of  the  condition  code 
values  and  the  meaning  of  each: 


CONDITION 

CODE 

INDEX  Error  Message 

1  The  starting  point,  ST,  is  outside  the 
input  data  array.  (Too  large) 

2  The  starting  point,  ST,  is  outside  the 
input  data  array.  (Too  small) 

3  The  input  data  variables,  ST,  RNG,  are 
improperly  defined  and  result  in  a 
condition  which  places  the  stopping 
point  before  the  starting  point. 

U  The  user  supplied  value  of  TH  is  less 

than  TL  when  NAUTO  =  1 . 

b.  The  subroutine  PACK  is  uded  to  convert  the  fixed 
point  output  data  into  a  form  suitable  for  trans¬ 
mission  to  remote  terminal  via  TSS. 


c.  If  an  output  file  code  (IFCODE)  is  not  specified 
then  no  data  will  be  transferred  and  control  will 
immediately  be  transferred  back  to  the  calling  pro¬ 
gram. 

d.  Flow  Chart:  Page  9-82 

e.  Cross  Reference  Table:  Page  9-218 
THEORY  OF  OPERATION 

a.  Upon  entering  the  subroutine  the  input  parameters 
are  first  tested  to  verify  that  the  range  of  the 
desired  plot  is  within  the  limits  of  the  input  data 
array,  DV. 

b.  If  NAUTO  =  0,  the  input  data  is  scanned  to  determine 
the  maximum,  TH,  and  minimum,  TL,  values  to  be 
plotted.  From  these  two  variables  the  upper  and 
lower  limits  of  the  plot  are  determined. 

c.  For  magnitude  data,  TH  and  TL  are  tested  to  deter¬ 
mine  if  the  input  data  is  bipolar  or  unipolar. 

d.  From  TH  and  TL  the  dependent  variable  range  of  the 
plot  is  computed  and  the  LSB  and  BIAS  to  be  used  in 
the  conversion  from  floating  point  to  fixed  point 
is  computed. 

e.  The  following  plot  parameters  are  processed  as  one 
record : 

1.  TH  -  Maximum  value  of  input  array 

2 .  TL  -  Minimum  value  of  input  array 

3.  IVST  -  First  value  of  input  array  which  was 

processed  for  output 

4.  NSTOP  -  Last  value  of  input  array  which  was 

processed  for  output 

5 .  OR  -  Independent  variable  value  for  first 

output  point 

6.  DEL  -  Independent  variable  spacing 

7.  RNG  -  Independent  variable  range  of  the 

output  data 

f.  The  following  plot  parameters  are  process*' a  for 
output  as  one  record: 
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1 .  BIAS 

2.  XLSB 

3. 

U.  ith 

5.  ITL 

6.  NOU'f 


Number  subtracted  from  input  data  before 
conversion  to  fixed  point 
LSB  of  output  data 
Not  used 

Largest  dependent  variable  output  value 
Smallest  dependent  variable  output  value 
Number  of  output  samples 


g.  The  dependent  variable  values  are  processed  for 
output.  The  mechanization  equation  for  this 
conversion  is  the  following: 

IDAT=(DV(J )-BIAS)/XLSB  IVST <  J  £ NSTOP 

The  integer  IBAT  is  then  processed  by  the  routine 
PACK  which  converts  it  into  2  characters  and  packs 
it  into  the  array  line  which  represents  our  output 
record  and  has  the  following  structure: 


3U  dependent  variable  samples 


header  character  pair 


trailer  character  pair 


6-31 


4^6  O 


bubfcUUIlNt  H  L I  b  T  <  C' V  I 


4.4  t>  l 
220L 

4  ?6* 

*.*:  #u 
*2/1 

22  /i 
22Ts 
**  H 
2.  /o 
22  r# 

22  K 
22  1* 
u'a 

tJ  wl 
22  t/ 

4  2~t 


ClMMuh/bL*  i/  bM  t  /  uC  I  »V  It  f  *0>  *  Via  I  bt*  >  •  ltt>l  I  tu  I 


t  wUA  VALlNC  l 

1 

i-  k  n 

b*  I  • 

b  T 

:HII 

<  u  1 

• 

1 

(Rll 

f  ?  )  . 

Lf 

1  .  1  HUI 

o_.  1 

• 

1 

*  Ml 

^  ‘  •  I  t> 

NAtJl  (, 

l.(  IM( 

• 

1 

.-'Ml 

06  5  » 

T  L 

1  .  1  ■  1  1  1  1  J> 

KMu 

NbK  ^ 
lM 

IKl.Ul 


lilMtllMlh  ('VI  1  I  »t  1,'a  I  It'  i 


uM  A  1  V&/-3  -  A  ,m/  *  .«/  i  „ui  0  2 


L'A  I  A  4.  b  f  *..  M/l  Vo  1  ?«  UOOUlyM/ 


If  t  ItiMut  •L-i-l-J  t»uTl>  tCO 
u  i  »^klU  (  i*  *  y  - 1 
Mil.  =  im.UlA'C'lv.tU 
Lf  -  IVIMS'.I 
vl  L  -  l  V  4  f«  i  **.  I 


I  1  J  .  I  Nl» 


I  -  C  \<  ki l  LMI  Li. A  I  I  Ml  1L  -  i  I 
k  t  *  ...  I  ..  *»<  «vi  f  I  i  f  T  l:  *'K/ 
hNi.-f.iii.-i.Kt:,  i 


•.i  K  ft.M  xl A\ 

It  (  L.u  ,IJ  I  v.,  H  7*. 

I-  M*  -  I  .  -  ‘  I 


(  ...  I  i  Mi  l 

i . i  1 1  » >- :  r  :  -I  »  >  ■  5  •  :  ^  i 
>.i  i •  i 


22*i 

It  (  1  V3I.OC  .f»ni  )  Nf-  ^  t  L  -  1 

22«< 

NSU.f  -  lHX««ST-L'R*RNCI/OtLI 

22*-> 

tk  =  M-*UVS1-l>*utL 

22  **♦ 

I  h  (  NS  lor1  .1/  1  .N1 1  0  )  NSTOP  -  N1TL 

22*t> 

IMNiTL+.Lt.JVST)  NCrPl  =  - 

22*o 

1 1- (NtTM-.Nt- .  u)  0>C‘  TO  1101/ 

22W 

lKlt.Ll.il  U  TO  4u 

22*- 

U'_  •»*  J-  1  Vb.  T  .NS  TLI- 

22** 

AsAMAAilAtSloVIJl  •  .XMlNOfc.  ) 

2300 

It  ILt  .uO.c  )  i.O  TL  4l 

23o  1 

L»V  l  J  1 1 1  0.0*  At  UO  lo<  x  • 

230*. 

ot,  To  a: 

230* 

.1  0Viol~Ai.lL/10IX) 

2  <  On 

./  LONUNOc 

<130  5' 

4  S  OCNTiM'A 

2->oo 

->u  COM  JNLtt 

i.iu? 

in  =  t/Vi  i  »sn 

2301 

11  =  lfi 

230* 

It-  INC  It  C  .M  .o  1  L/C  TL  lloo 

23  lo 

DO  11  It  J-  J  VS  I  .NSTCP 

2311 

i,t  c  ap.AX  It  oV  1  J  1  .  To  1 

23l2 

It  -  AKlNliLtVIJl.Tcl 

1<  (■  UM  iMJl 

23  i* 

it  (In.ol  .TL  )  (>t  TO  ?<t01 

til*. 

NCIt.O  =  S 

23  lo 

L.O  Tt  lloo 

2317 

ivUl  COM  INC t 

231*. 

otL  =  i/tO*l  cl  AT  INSXT’1 

2->l* 

NUlTHNSH  P-IVST1/NSK.P*! 

ufc/  I  1/  7  j 


iNf'Ut  Ll^IiNu 


AlIuFLLm  CHAKl  ill  -  t-ML/^CL  KAUilM 


CAKD  Nb  •••• 


LINUNTS 


4*dl- 

zy<  i 

<£3.\> 

*.’A.  '> 

4.  jjI 

4  *3* 

i  jj*» 

«.  i.W 
«\*3- 

/-»*»! 
*  .w. 

..  .»•*  . 

»'  '•♦*♦ 


L  INI.  Cl  1-btU  I  Ihl 
1 1N».  t^l-bl'UL  <  U  I 

LINlI 31-lVil 
L1nl(*i  UNiltll 
lINi I b l=olLL (UK ) 

LiNt  (  oUcLU  ll>t L  1 
LlNc  C  7  Hl.uUL  CkNO  I 
Ml  1  I L I  1! L  LUt  1  lINf 
1ML1.1U.  1.  ul.Ll.ll..  31  LOT'.  7t-u 
11 ClH.Lt.u.ul  1LON-1 
11  ( 1 L  .Ut  . u .U )  1C.LN-. 

HllM.bl.u.u.ANL.TL.lT.iuOl  ILLS-  j 

Mrls  AHtl  1  (*.o|  Tm  )  1*1  il  1L  I  J 

!lL-  ALuuluCllHl 

UlM  -  lllXCilLl 

*1  II  Il.l  I  .0.4-1  LU  1C  1!»0 

x-l  IH/I  10.1. ••!  it.;>i  1 

bt  It.  till' 

I'll.  X-l  Il»(  tl/.tiMUlUM  III  il  HI  II 

i(.L  UNI  1  NUl 

N  -  .  i  . 

II  CA.Ot.l.L.  tNtl  .X. Lb. 1.11  llst 

I I I  a  .  1. 1  •  t  *  v,  .  i*  l.u  .  X  .  1. 1  .  i  .  M  N  -  ** 
lllt.il.i.'  .  .iil  .X  .LI  .h.1  I  N  —  * 

11  (  K  .I.  I  . *♦  .  !. .  Al.l i .  X  .  I  l  .  7  . i*  I  N  —  l 
1  .  It.bl.  /  .I'.Al.L'.X.Lr  .W.O  IN=  1  l 
illilt.il.i  J  I  t,i  li.  _> ■  « > 

mi  tri.  <  i  ci.i*i  1 1  .i'i»ni  m  i 


£-33 

•  ••• 


i  l  jc  i 


****** 


.>  >  v 


Ml, AT  IHl  /lli'.l<«*tl*b»  « 


Cf-ia*- 


2JU' 

llM  lt.ut 

Zl'Zk 

»LbL  =  MH/40i>l-.U 

t  jt<- 

hi  tb=ll.O 

l,U  1  Ut  jOV.hoI  t4h2  t 

AMX  -  O.U 

AMN  --TH 

<! 

AMO  »'lh/t.U 

ci!>7 

ulISC-  Th/ioo.O 

Ou  10  *tO-4 

2^t?v 

<*u  1 

A MX  -  TH 

jou 

AMU  -  0.0 

AMli  -  lM/2.< 

«:  J04. 

I/INL-  1  h/  i  l  l  .O 

AU  . 

Ol  lo  IKI 

J  :><•*♦ 

H' 1 1. 

AM.X  -  V’. 

je>> 

AMU  =•-!> 

<?JCC 

AM.  -  (• 

4.  jt  / 

uiNC-  ll/OO.u 

*  .tot 

l;y  lo 

i.'  7v 

*  ;>  o 

<_  i  7** 
/-■»  /  > 
4'^  ^ 

..  s  f  i 


tv  L.  MML-/.bbUf"U) 

M  -  (  i  f  i  X  Um'/U'ibM*^ 
iY  <*Mu«v#1  •  iv»o.u  I  N-loO 

MtvX  *  (MUTin  .(> 

lh  UKA.LT.  Ini  *MX  =  *MXMO.L‘ 

Aftf* -  ^M>-rLUUN| 

ir  I  •  I  l  •  Aul»«N  »LT  •  10C>  I  li=N*  U> 

i  f  IN.t  w  •  *  o  I 

4 1  * ' 


1 


2>i  t. 

iJwU 

i 

4  3  jC 
tJl  / 

ZKt 
/  Jv  >* 

<.  ,>**i» 

<;.^4 

4._>V> 

O^- 

<Mv7 

*:-»w 

«^'>v 

4  ^t|J 

4.*tO< 

<_««  w 

f|l  •» 


lHN.tW.t>i>.liK.h.tU.7«»l  N*LO 
tlAi=  aNX'HuAIINI 

XLib-li.l'i' 

IMN.Lc.^ul  XLi8=0»i*l 
lMN.i.I.'.O.ANG.N.t-t.boJ  XLSB=b.U* 

-,t,i  lHNirkL.Nt.llol  To  11W 
HL-(lL-l.iAS)/XLbt 
lTh=-tU.-blA8J/XL8L 
Llljt  (  1  I  i-bl'OLt  i'-I  AS  I 
Lint  ( <•  >  -  ■>!  UlIXLSOI 

L  i  N t  (jUK 

Lltit  «-.!=■  i  Ih 

Lilli  (  k i  =  i Tl 

LlULI  i*  I=i«l0i 

Kk  i  tl  (  Hill'.-  >  LlNi. 

Lit.l  (  1  >-lSUH 
c  Jilt  tic  >-LLlCh 
oi.;  iwt-i 

luil=]t 

mil  kmCX  41-ijiT  flWUtlPlTtLlUttllliiOl 

|,L  <.1,1  J^iVSI.NMlk.NSKk 

li,  A1  -(cVI  i  l-t  iAil/XLSM 

i  1 1  li  A 1  »L  i  *-*.li**k  )  1  uAT  =—• *livk 
LALL  HLMJl.Al,IW1.tIb!T,Llr<L,»i<..» 

V  L  II  M'U 

i  1. 1  nlilLl.KLUl  LlNt 

1*L'-  i 
I«.i  l^l« 


1 


SUBROUTINE  RTOPDB 


1  .  MODULE  IDEBTIFICATION: 


Name 

Classification  Code 

Reference  Number 

RTOPDB 

Peripheral 

103 

RTOPM 

Peripheral 

110 

RT0PM2 

Peripheral 

111 

XYTOM 

Peripheral 

105 

XYT0M2 

Peripheral 

106 

XYTODB 

Peripheral 

10U.108 

2 .  PURPOSE: 


This  subroutine  performs  a  rectangular  to  polar 
conversation  on  the  input  arrays.  Depending  on  the  entry 
point  used,  the  output  is  in  two  arrays  in  the  form  of 
linear  magnitude  and  phase  (RTOPM),  magnitude  squared 
and  phase  (RT0PM2)  or  magnitude  in  decibles  and  phase 
(RTOPDB);  or  in  one  array  of  linear  magnitude  (XYTOM), 
linear  magnitude  squared  (XYT0M2)  or  magnitude  in  de¬ 
cibels  (XTTODB). 


3.  INPUT  PARAMETERS:  None 

U.  CALLING  SEQUENCES: 

Phase  and  decibel  magnitude 


CALL  RTOPDB  (X,Y,M,P) 


Where:  X 

Y 

M 

P 

M(J)  =  10.0  *  ALOGIO 
P(J)  =  ~8^— °  *  ATAN2 


contains  the  Input  Waveform  -  R 
contains  the  Input  Waveform  -  I 
contains  the  Output  Waveform 
contains  the  Output  Waveform  -  P 
[  X(J)  *  X(J)  4-  Y(J)  *  Y(J)  ] 
(Y(J),  X(J)) 


Phase  and  linear  magnitude 
CALL  RTCPM  (X,Y,M,P) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
M  contains  the  Output  Waveform  -  M 
P  contains  the  Output  Waveform  -  P 


M(J)  =  V  X<J)  *  X<J)  +  Y<J>  *  Y<J> 

P(j)  «  180^0  *  ATAN2  (Y(J),X(J)) 

TT 

Phase  and  magnitude  squared 
CALL  RT0PM2  (X,Y,M,P) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
M  contains  the  Output  Waveform  -  M 
P  contains  the  Output  Waveform  -  P 

M(J)  =  X(J)  *  X(J)  +  Y(J)  *  Y(J) 

P(J)  «  *  ATAN2  (Y(J),  X(J)) 

Linear  magnitude  and  no  phase 
CALL  XYTOM  (X,  Y,  M) 

Where:  X  contains  the  Input  Waveform  -  R 

Y  contains  the  Input  Waveform  -  I 
M  contains  the  Output  Waveform 

M(J)  =  yjX(J)  *  X(J)  +  Y(J)  *  Y(J) 
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Magnitude  squared  and  no  phase 
CALL  XYT0M2  (X,Y,rf) 


Where :  X  contains  the  Input  Waveform  -  B 

Y  contains  the  Input  Waveform  -  I 
M  contains  the  Output  Waveform 

M(J)  =  X(J)  *  X(J)  +  Y(J)  *  Y(J) 

Magnitude  in  dB  and  no  phase 

CALL  XYTODB  (X,Y,M) 

Where:  X  contains  the  Input  Waveform  -  B 

Y  contains  the  Input  Waveform  -  I 
M  contains  the  Output  Waveform 

M(J)  =  10.0  *  ALOGIO  [x(j)  *  X(J)  +  Y(j)  *  Y(j)] 

RESTRICTIONS.  BgggIBgMjjMTS .  MISCELLANEOUS  LATA 

a.  For  entry  points  BTOPDB  and  XYTODB  the  minimum  out¬ 
put  value  is  -200  dB. 

b.  The  range  of  the  output  phase  angle  is  from  -190°  to 
+170°.  The  break  point  where  phase  "rolls  over"  was 
chosen  to  be  170°  instead  of  180°  because  many  yave- 
foxms  have  a  large  number  of  samples  around  180°  and 
each  time  the  waveform's  phase  changes  from  180  -B 
to  180  +E  the  output  phase  would  change  from  idO-B 
to  -180+E.  This  results  in  a  phase  plot  with  many 
undesirable  and  unnecessary  discontinuities. 

c.  Flow  Chart:  Page  9-156 

d.  Cross  Reference  Table:  Page  9-228. 
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SUBROUTINE  SCANNR 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 


SCANNR  PERIP  313 

2.  PURPOSE: 

This  subroutine  scans  an  antenna  pattern  in  the  form 
of  an  input  array  and  determines  the  location  and  gain  of 
(1)  the  main  lobe,  (2)  the  highest  ten  sidelobes,  and  (3) 
the  first  twenty  sidelobes  on  either  side  of  the  main  lobe. 
In  addition,  the  first  minimum  on  either  side  of  the  main 
lobe  is  determined. 

3.  INPUT  PARAMETERS: 

None 

4.  CALLING  SEQUENCES: 

CALL  SCANNR  (DV) 

Where:  DV  contains  the  Input  Waveform  in 

dB  gain 

5.  RESTRICTION,  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  All  output  data  from  this  subroutine  is  in  the 
form  of  a  printed  listing  which  is  formatted  for 
ease  of  interpretation. 

b.  If  a  vestigial  lobe  exists,  the  minimum  between  it 
and  the  main  lobe  will  be  detected  as  the  main 
lobe  null. 

c.  Flow  Chart:  Page  9-161 

d.  Cross  Reference  Table:  Page  9-229 
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6. 


THEORY  OF  OPERATION 


a.  The  antenna  pattern  data  is  scanned  from  left  to 
right  (increasing  0)  and  the  gain  maxima  are  stored 
in  a  22  element  array.  Within  this  array,  the 
maxima  are  ordered  according  to  gain  with  the  high¬ 
est  in  element  1  and  the  lowest  in  element  22. 

b.  The  pattern  data  to  the  right  of  the  main  lobe  is 
scanned  to  locate  the  main  lobe  null  and  the  first 
twenty  sidelobes.  After  this  data  is  processed 
and  printed,  a  third  scan  is  made  to  the  left  of 
the  main  lobe  for  the  first  null  and  the  first 
twenty  sidelobes. 

c.  The  logical  equation  used  to  determine  the  peaks 
of  the  lobes  is  as  follows: 

Given  three  adjacent  points  nV(NTHPT-l),  DV(NTHPT), 
and  DV(NTHPT+1)  the  following  difference  parameters 
are  computed. 

DELDV1  =  DV(NTHPT)  -  DV(NTHPT-l) 

DELDV2  =  DV(NTHPT+1)  -  DV(NTHPT) 

then  if  DELDV1  >  0  and  DELDV2  <  0, 

the  point  DV(NTHPT)  is  considered  to  be  a 

maxima,  i.e.,  the  peak  of  a  lobe. 

d.  In  order  to  find  the  nulls  of  the  main  lobe,  the 
difference  parameters  defined  in  6(c)  are  tested 
for  the  following  condition: 

DELDV1  <  0  and  DELDV2  >  0 

If  this  condition  is  satisfied  then  DV(NTHPT)  is 
considered  to  be  a  null.  After  the  first  null  is 
found,  a  flag  is  set  which  causes  the  null  test 
to  be  bypassed  for  the  remainder  of  the  scan  on 
that  side  of  the  main  lobe.  It  should  be  noted 
that  the  null  logic  will  consider  the  minimum  be¬ 
tween  the  tpain  lobe  and  a  vestigial  sidelobe  (if 
one  exists)  to  be  the  main  lobe  null. 
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SECTION  7 


SUPERVISORY  MODULES 


This  section  includes  those  modules  that  supervise 
existing  stimulus/transfer  function  modules  to  simulate 
the  frequency  scanned  and  time  scanned  array  radar  systems. 
This  group  includes  ANTARY  and  TSRPAT,  which  are  located 
in  Volume  III. 
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SECTION  8 


SUBORDINATE  MODULES 


This  section  includes  all  those  modules  which  are  sub¬ 
ordinate  to  the  larger  stimulus /transfer  function  modules, 
peripheral  modules  or  supervisory  modules.  They  include 
both  subroutines  and  functions.  Included  in  this  group  are 
the  following  modules: 


ABORT 

IBOOL 

ANTINT 

IFLD 

AZGA IN/ELGA IN 

IPACK 

BLOCK  DATA 

PACK 

DBLKX 

RRAND 

The  following  subordinate  modules  are  related  to  the 
bistatic  target  model  and  target  imaging  and  are  located 
in  Volume  I,  Part  3: 

BISTGT 

GAM 

BESS 

EXPT 

EXPI  is  also  located  in  Volume  IV,  Part  2. 
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SUBROUTINE  ABORT 


1 .  MODULE  IDENTIFICATION: 

Name  Classification  Reference  Number 

ABORT  Subordinate  Not  user  referenced 

ERRMSG  Subordinate  Not  user  referenced 

2.  PURPOSE: 

This  subroutine  is  used  to  print  coded  error  messages 
to  the  user. 

3 .  INPUT  PARAMETERS: 

Name  O/R  T  Description 

NCODE  R  I  Error  message  code  number. 

0  <  NCODE  <  999 

4 .  CALLING  SEQUENCES: 

CALL  ABORT  (NCODE) 

A  message  is  printed  which  states  a  fatal  error  has 
occurred  and  execution  will  terminate. 

CALL  ERRMSG  (NCODE) 

A  message  is  printed  which  states  a  non-fatal  error 
has  occurred,  a  fix-up  procedure  done,  and  execution 
continued . 

5.  RESTRICTIONS,  REQUIREMENTS,  MISCELLANEOUS  DATA 
a.  Flow  Chart:  Page  9-81 

h.  Cross  Reference  Table:  Page  9-218 
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SUBROUTINE  ANTI  NT 


1 .  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

ANTINT  Subordinate  (301) 

2.  PURPOSE: 

This  subroutine  generates  an  interpolation  table  from 
the  user  provided  sampled  antenna  pattern  data.  The  inter¬ 
polation  table  is  subsequently  used  by  the  functions  AZGAIN 
and  ELGAIN  to  compute  antenna  gain. 

3.  INPUT  PARAMETERS : 

Name  0/R  T 

NPT  R  I 

BSIT  0  F 

ANTP( 2,75)  R  F 

4.  CALLING  SEQUENCES: 

CALL  ANTINT  (NPT,  BSIT,  ANTP ,  CDEF) 

Where:  ANTP  contains  the  Input  Antenna  Pattern 

COEF  contains  the  Output  Interpolation 
coefficients  table. 

5.  RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  If  no  antenna  pattern  is  specified  then  the  co¬ 
efficients  table  is  set  up  to  represent  an  omni¬ 
directional  antenna  of  unity  gain. 


Description 

Number  of  user  specified  antenna 
gain  values 

Angle  between  zero  angle  of 
user  specified  pattern  and  zero 
angle  in  radar  coordinates 
(usually  =  0.0) 

User  specified  gain  pattern  array. 
ANTP(1,J)  is  the  gain  of  the  Jth 
sample . 

ANTP(2,J)  is  the  angle  of  the 
Jth  sample. 
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b.  Flow  Chart:  Page  9-140 

c.  Cross  Reference  Table:  Page  9-226. 
THEORY  OF  OPERATION 


6. 


This  module  implements  a  simplified  antenna  represen¬ 
tation  which  is  based  on  computing  the  antenna  gain  as  the 
product  of  samples  from  an  azimuth  and  elevation  pattern 
cuts  through  the  peak  of  the  main  lobe.  The  azimuth  pattern 
is  true  antenna  gain  in  power  units.  The  elevation  pattern 
is  normalized  to  unity  gain  at  the  peak  of  the  main  lobe. 

The  parabolic  interpolation  scheme  mechanized  herein 
is  described  in  the  following  paragraph.  Given  three  points 
a  parabolic  curve  can  be  determined  which  passes  through 
them. 


where  A,  B,  and  C  are  the  coefficients  to  be  determined 


XYl  2  -  XY13 
X2  -  X3 


B  =  XY12  -  h  *  (XI  +  X2) 

C  =  Y1  -  (X1A  +  B)  *  XI 

Y1  -  Y2  Yl  -  Y3 

XYl 2  =  xi  '-"X2  and  XY13  “  XI  -  X3 

In  determining  the  coefficient  table  elements  a  double 
parabolic  coefficient  is  used  where  possible.  This  is  done 
by  determining  2  sets  of  coefficients  for  2  curves  through 
4  points  and  averaging  the  coefficients.  Typically  a  double 
parabolic  curve  fit  is  smoother  than  a  single  parabolic 
curve  fit. 
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SUBROUTINE  AZGAIN/ELGAIN 


1.  MODULE  IDENTIFICATION: 


Name 

AZGAIN 

ELGAIN 

2.  PURPOSE: 


Classification  Code 
Subordinate 
Subordinate 


Reference  Number 
Not  user  referenced 
Not  user  referenced 


This  function  is  to  compute  the  antenna  azimuth  or 
elevation  gain  for  any  angle  within  the  range  of  the  user 
specified  antenna  data. 

3.  INPUT  PARAMETERS: 

Name  O/R  T  Description 

ANTANG  R  F  The  angle  of  the  antenna  main 

lobe  in  radar  coordinates 

4.  CALLING  SEQUENCES: 

CALL  G  =  AZGAIN  (AZANGL)  *  ELGAIN  (ELANGL) 

Where:  GA  is  the  antenna  gain  corresponding  to 

an  azimuth  angle  of  AZANGL  and  eleva¬ 
tion  angle  of  ELANGL. 

5.  RESTRICTIONS.  REQUIREMENTS  AND  MISCELLANEOUS  DATA 

a.  The  coefficients  table  must  have  been  successfully 
initialized  by  the  subroutine  ANTINT  before  the 
function  can  be  executed. 

b.  Flow  Charts:  Pages  9-159/160 

c.  Cross  Reference  Tables:  Page  9-229. 

6.  THEORY  OF  OPERATION 


The  input  angle  ANGL  is  tested  to  determine  if  it  is 
within  the  range  of  the  input  data.  If  it  is  outside  the 
range,  the  gain  value  returned  is  zero.  Otherwise,  the 
gain  is  computed  from  the  gain  coefficients  (A,B,  and  C) 
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that  are  closest  In  angle  to  ANGL. 

AZGAIN  =  (ANGR  **  2)  *  A  +  ANGR  *  B  +  C 

Note:  AZANGL  is" the  angle  measured  in  degrees  CCW  from 
any  desired  reference  direction.  ELANGL  is  the 
angle  measured  in  degrees  up  from  the  horizon. 
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SUBROUTINE  BLOCK  DATA 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

BLOCK  DATA  Subordinate  None 

2.  PURPOSE;  ' 

This  non-executable  module  initializes  the  labled 
common  BLKRND  which  contains  the  random  number  array  used 
by  the  random  number  generator. 

3.  INPUT  PARAMETERS: 

None 

4.  CALLING  SEQUENCES: 

None 

This  module  is  used  only  at  load  time  to  initialize 
the  labeled  area. 

5.  RESTRICTIONS.  REQUIREMENTS,  MISCELLANEOUS  DATA 

a.  The  numbers  used  to  initialize  the  random  number 
array  were  derived  from  tables  of  random  numbers 
published  by  the  RAND  Corporation:  A  Million 
Random  Digits  with  100.000  Normal  Deviates. 

the  Free  Press  of  Glencoe,  New  York,  1955. 

b.  Flow  Chart:  Page  9-124 

c.  Cross  Reference  Table:  Page  9-224 
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SUBROUTINE  DBLKX 


1.  MODULE  IDENTIFICATION; 

Classification  Code  Reference  Number 
Subordinate  Not  user  referenced 

Subordinate  Not  user  referenced 

This  subroutine  is  used  to  transfer  a  defined  block 
of  data  from  one  array  to  another  array. 

3.  INPUT  PARAMETERS: 

a .  DBLKX 

Name  0/R  T  Description 

NST  R  I  First  element  of  the  input  array 

to  be  transferred.  Also,  the 
first  element  of  the  output  array 
to  be  used. 


Name 

DBLKX 

SDBLKX 

2.  PURPOSE: 


NWORD  R  I  Number  of  elements  to  be  transferred. 

b .  SDBLKX 

Name  0/R  T  Description 

NST  R  I  First  element  of  the  input  array 

to  be  transferred. 

NWORD  R  I  Number  of  elements  to  be 

transferred. 

LOC  R  I  First  element  of  the  output  array 

to  be  used. 

4.  CALLING  SEQUENCES: 


CALL  DBLKX  (NST,  NWORD,  X,  Y) 

Where  X  contains  the  Input  Waveform 

Y  contains  the  Output  Waveform 


Y(NST+J)  -  X(NST+J)  J-l,  NWORD 
CALL  SDBLKX  (LOC,  NST,  NWORD,  X,Y) 


Where:  X  contains  the  Input  Waveform 

Y  contains  the  Output  Waveform 

Y(LOC+J)  -  X(NST+J)  J-l,  NWORD 
5.  RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  All  input  data  to  this  subroutine  is  provided 
through  the  calling  argument  list.  This  sub¬ 
routine  is  not  directly  available  to  the  program 
user,  but  is  used  by  other  subprograms  for  data 
transfer. 

b.  Flow  Chart:  Page  9-122 

c.  Cross  Reference  Table:  Page  9-223. 
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FUNCTION  I BOOL 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Name 

IBOOL  N/A 

2.  PURPOSE: 

This  function  was  used  when  checking  out  RADSIM  on 

IBM  computers.  It  has  been  replaced  by  the  standard 
Honeywell  intrinsic  function  B001. 

3.  INPUT  PARAMETERS:  N/A 

4.  CALLING  SEQUENCES:  N/A 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-119 

b.  Cross  Reference  Table:  Page  9-223 
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FUNCTION  IFLD 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

IFLD  N/A 

2.  PURPOSE: 

This  function  was  used  when  checking  out  RADSIM  on 
IBM  computers.  It  has  been  replaced  by  the  standard 
Honeywell  intrinsic  function  FLD. 

3.  INPUT  PARAMETERS:  N/A 

4.  CALLING  SEQUENCES:  N/A 

5.  RESTRICTIONS.  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  Flow  Chart:  Page  9-123 

b.  Cross  Reference  Table:  Page  9-224 
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FUNCTION  I PACK 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

IPACK  Subordinate  Not  User  Referenced 

2.  PURPOSE: 

This  function  subprogram  is  used  to  pack  two  or  more 
groups  of  bits  into  a  standard  computer  word  (36  bits). 

3.  INPUT  PARAMETERS: 


Name 

Q/R 

T 

Description 

ISTR 

R 

I 

The  lowest  bit  into  which  data 
transfer  is  to  occur.  The  bit 
number  is  counted  starting  with 
the  LSB. 

ISTR  =  0  corresponds  to  the  LSB  and 
ISTR  =35  corresponds  to  the  MSB. 

IDATA 

R 

I 

The  data  which  is  to  be  packed 
into  the  output  word. 

IWORD 

R 

I 

The  word  into  which  the  data  is  to 
be  packed. 

4.  CALLING  SEQUENCES: 

CALL  IOUT  -  IPACK  (ISTR,  IDATA,  IWORD) 

I  OUT  -  IWORD  ©  ( IDATA*  (2**ISTR))  ©  -  inclusive  or 

5.  RESTRICTIONS  t  REQUIREMENTS.  MISCELLANEOUS  DATA 

a.  The  field  (bit  locations)  of  the  output  word  which 
is  to  receive  the  input  data  must  contain  zeros 
before  execution  of  the  function. 

b.  If  the  data  to  be  placed  in  the  output  word  is 
larger  than  the  allocated  field  then  either  the 
data  in  a  higher  field  will  be  destroyed  or  if 
it  is  the  highest  field  then  some  of  the  most 
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significant  bits  of  IDATA  will  extend  past  the 
standard  word  length  and  be  lost. 


c.  Flow  Chart:  Page  9-118 

d.  Cross  Reference  Table:  Page  9-223 
6.  THEORY  OF  OPERATION 

The  input  data,  IDATA,  is  multiplied  by  the  power  of 
two  equivalent  to  the  number  of  bit  positions  that  it  must 
be  shifted  in  order  to  line  up  with  its  field  in  the  output 
word.  Once  the  data  has  been  shifted  to  the  proper  posi¬ 
tion,  then  it  is  merged  into  the  output  word,  IWORD,  by 
using  an  inclusive  OR  operation. 

This  function  is  not  directly  available  to  the  simu¬ 
lation  user,  but  is  used  by  subprogram  within  the  simulation. 
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SUBROUTINE  PACK 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

PACK  SUBORDINATE  Not  user  referenced 

2.  PURPOSE: 

This  subroutine  takes  a  word  in  integer  format  and 
converts  the  lower  13  bits  of  the  word  into  two  8  bit  ASCII 
characters.  The  ASCII  characters  are  subsequently  trans¬ 
mitted  to  a  remote  plotting  terminal  via  TSS.  Remote  plot¬ 
ter  accuracy  assumed  to  be  12  bits  maximum. 

3.  INPUT  PARAMETERS: 


Name 

0/R 

T 

Description 

IWD 

R 

I 

Output  array  word  location  pointer. 

IBIT 

R 

I 

Bit  displacement  to  be  used  in 
placing  the  characters  in  the 
output  word 

4.  CALLING  SEQUENCES: 

CALL  PACK  (IDAT,  IWD,  IBIT,  IARY,  $1111) 

Where:  IDAT  is  the  Input  word 

IARY  is  the  Output  array 

IIII  is  the  statement  number  to  which  con¬ 
trol  is  transferred  if  the  output 
array  is  full  (68  characters) 

5.  RESTRICTIONS,  RECOMMENDATIONS.  MISCELLANEOUS  DATA 

a.  The  maximum  array  length  is  restricted  to  68 
characters  to  be  compatible  with  TSS. 

b.  The  maximum  allowable  accuracy  of  the  word  is 
13  bits,  i.e.  bits  0  through  24  must  be  zero. 

c.  Flow  Chart:  Page  9-88 

d.  Cross  Reference  Table:  Page  9-219 
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THEORY  OF  OPERATION 


The  input  word  IDAT  is  shown  pictorially  as  follows: 


IDAT  o 


Bit  0  is  the  MSB  and  Bit  35  is  the  LSB. 


This  word  is  divided  into  two  characters,  II  and  I2as 
illustrated  by  the  following: 


1 

l 

I 

l 

s 

s 

1 

3 

l 

i 

3 

3 

3 

3 

II 


12 


Next  12  is  tested  to  determine  if  any  TSS  command  charac¬ 
ters  have  occurred.  If  so,  the  value  of  12  is  incremented 
by  +1  or  -1  as  follows: 


TSS  Character 

Actual  No. 

Increment 

DC1  (XON) 

218 

+1 

CAN 

00 

o 

ro 

41 

EOT 

48 

41 

RUBOUT 

177g 

-1 

SOH 

1 

-1 

ETX 

3 

-1 

Since  II  is  6  bits  in  length  a  rubout  (177g)  cannot 
occur.  Therefore,  the  remaining  TSS  characters  to  protect 
against  have  values  which  are  less  than  32-jp*  Accordingly, 
II  is  tested  and  if  its  value  is  less  than  32^0  then  bit  7 
of  II  is  set  "on" to  preclude  occurrence  of  the  undesired 
character . 

Finally,  II  and  12  are  packed  into  the  output  word 
pointed  to  by  IWD.  The  parameter,  IBIT,  specifies  the  dis¬ 
placement  of  II  from  bit  0  of  the  output  word,  IARY  (IWD). 


8-19 


ilihkUll  lht  HACK  ( I  CAT  •  1  Mb  ,  1  bl  T  ,  1  AkV  .*  I 
C  A  I  A  IKCl/cl  /7/,lftlT7/Ul«(>/,  1XCN/CV1/,  ILANL/CJli/ 
UATa  1 1  U I /*./,  1  C/017?/,  1C  UN  A/  1/,  I  CGNC/i/ 
UIMlNCICN  lAKY(l) 


I- 1  Ak  Y  l  i  H  L)  ) 

1  1-H.OI  ZJ  1 1-,  ]  LAl  1 

1<:  =  h  ci<(  .'v,  /  ,  Ida  1 1 

lHli.ll  .i.-.ANi-.I  1  .GT.I.)  11-1  1*IHI  7 

II-  (  K-.Lu.  IK.1  .CK  .  1  /  .tU.lCONA.Uk.l?.LW.  lCciNc  I  U  =  lc~l 
lH|J^.tl..lXU.'»..JK.i:,.!A.lCANC.Ul.li.LC.ltUl>  1  ?-l i«l 

lMI->iwkl  1C  *  1 1  *  1 


Hl'l  1 1  -  a  ?  iVi  1AKY1  IWCI  1-1  1 


1  !  11  -  1  I  l  I  ♦*' 


HIM  11.1  I  .v  ,1  aKY  UWI  1  )-  1  «■ 


if.  I  1  =  1  hi  I  ** 

If  Ijl  ll.ll.jM  ki  HlkN 


I  /  .  I  ill  Ail-  11  H( 


A«UH(WCHAkI  it.  1  -  MVC/iCC  kADSJH 


£1  M(  Mi 


1  Wl  -  1  Mi  -*  i 


1  k  I  1  m  .1  I  •  I  I  I  k  i  1 i.k  N  1 


li  I  ck  K 


FUNCTION  RRAND 


1.  MODULE  IDENTIFICATION: 

Name  Classification  Code  Reference  Number 

RRAND  Subordinate  101,  102 

(Parameter  initia¬ 
lization  only) 


2.  PURPOSE: 

This  function  generates  random  numbers  for  use  in  var¬ 
ious  subprograms  of  the  simulation  model.  Samples  from  the 
uniform,  Gaussian  and  Rayleigh  distributions  can  be  generated. 
In  addition,  two  statistical  target  models  are  incorporated 
which  are  based  on  the  Chi-square  distribution  with  2  and  4 
degrees  of  freedom. 

3.  INPUT  PARAMETERS:  (User  supplied  data  through  namelist 

NL101) 


Name 

O/R 

T 

Description 

MADl 

0 

I 

The  starting  address  for  select¬ 
ing  random  numbers  from  the 
random  number  array  (14MM>4128) 

UMEAN 

0 

F 

Mean  value  of  the  uniform 
distribution 

UUEXT 

0 

F 

Width  of  the  uniform 
distribution 

XMEAN 

0 

F 

Mean  value  of  the  Guassian 
dis  tribution 

SIGMA 

0 

F 

Standard  deviation  of  the  Guassian 
and  Rayleigh  distributions,  and 
average  cross  section  for  target 
models 

NTYPE 

R 

I 

Control  integer  which  specifies 
the  type  of  distribution  to  be 
generated 
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Name  0/R  T 


Description 

NTYPE  R  I  NTYPE  =  1  Uniform  distribution 

(floating  point  output) 
=  2  Rayleigh 
=  3  Gaussian 

=  4  Uniform  distribution 
(integer  0  to  2^-*'> 

=  5  Not  used 
=  6  Swerling  Target 
Models  #1  and  #2 
=  7  Swerling  Target 
Models  #3  and  #4 
=  8  Sine 

NRAND  R  I  Array  containing  random  numbers 

which  are  used  to  generate  the 
output  variates 

Variables  initialized  during  execution  of  the  initia¬ 
lizer  load  module  based  on  user  supplied  data.  These 
intermediate  values  are  generated  in  order  to  minimize 
execution  time  of  the  function,  i.e.,  to  avoid  re¬ 
peating  the  same  calculations  for  each  entry  to  RRAND. 

SIG2SQ  =2.0  * SIGMA  *  SIGMA 
UL  =  UMEAN  -  0.5  *  UUSXT 
UEXT  =  UUSXT/2  **  35 

CALLING  SEQUENCE : 

VAR  =  RRAND  (NTYPE) 

VAR  will  contain  the  random  sample  generated  by  the 
function  from  the  NTYPE  probability  distribution. 

RESTRICTIONS,  REQUIREMENTS.  MISCELLANEOUS  DATA: 

a.  Before  any  call  can  be  made  to  the  function  RRAND  th 
the  input  data  must  be  loaded  using  namelist  101. 

Any  subsequent  changes  in  random  variate  distribu¬ 
tions  are  also  made  through  namelist  101. 

b.  For  convenience  and  to  minimize  program  steps  the 
array  IRAND  was  equivalenced  to  the  array  NRAND 
but  displaced  by  one  location.  This  structure 
allows  an  address  of  zero  to  be  used,  i.e.,  an 


address  of  zero  will  access  IRAND(O)  which  overlays 
NRAND(1n.  If  this  were  not  done,  a  test  would 
have  to  be  performed  on  MAD1  to  ensure  that  an 
address  of  zero  did  not  occur. 

c.  A  list  of  random  numbers  suitable  for  initializing 
the  array  NRAND  are  included  at  the  end  of  this 
section. 

d.  Flow  Chart  Page  9-120 

e.  Cross  Reference  Table  Page  9-223 
6.  THEORY  OF  OPERATION: 

For  each  call  to  the  function  RRAND  a  number  IRND 
is  selected  from  the  random  number  table  IRAND.  The 
address  of  the  number  selected  from  the  table  is  MAD1 
which  is  also  a  random  number.  The  number  IRND  is  added 
to  the  random  number  JRND  which  was  generated  by  the 
previous  execution  of  the  function.  The  sign  bit  is 
set  to  zero  to  ensure  a  positive  number.  By  adding 
the  two  random  numbers  and  truncating  the  overflow,  a 
new  random  number  is  generated  which  is  also  called 
IRND.  IRND  is  placed  in  the  random  number  table  lo¬ 
cation  previously  occupied  by  the  original  IRND.  In 
this  manner  the  random  number  table  is  updated  by 
generating  new  random  numbers  and  inserting  them  in 
the  table.  From  this  random  number  IRND,  7  bits  are 
selected  to  determine  the  new  address  MADl  to  be  used 
in  the  next  call  to  the  function.  The  7-bit  address 
field  allows  the  addresses  to  range  from  0  to  127.  The 
random  number  IRND  just  generated  is  an  integer  having 
an  uniform  distribution  from  0  to  235  -  1.  Once  the 
random  number  is  generated  JRND  is  set  equal  to  IRND 
for  use  in  the  subsequent  executions  of  the  function. 

In  order  to  convert  this  number  to  a  floating  point 
number  r  having  a  uniform  distribution  from  0  to  1.0, 
the  following  conversion  is  used. 

r  -  FLOAT (IRND) / 2 ^ 5 

From  this  uniform  distribution  other  probability  dis¬ 
tributions  can  be  generated  by  using  transformations 
which  map  a  uniform  distribution  into  the  desired  dis¬ 
tribution.  The  following  is  a  list  of  the  transform¬ 
ations  used  in  the  function' 
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a . 

Uniform 

distribution 

a  ^  x  ^  b 

xn 

=  (b  -  a>rn  ■ 

+■  a 

/ 

b. 

Rayleigh  distribution 

P  (x )  =  *  J 

a 2 

xn 

=  In  rn 

1 

c . 

Gaussian 

distribution 

i. 

p(x^  =j2na 

(  x2\ 

■  W/ 


=  J 2 a2  In  r~  *  COS  2rrrn+1 


-(■£) 

d.Swerling  Target  Models  #1  &  #2  P(x)  *  _!  e  CT' 


xn  =  -a  In  rn 


o  =  average  cross  section 


This  distribution  charac terizes  a  target  with  a  large 
number  of  independent  scatterers  of  approximately 
equal  cross  section  (Swerling  cases  1  and  2). 

e.Swerling  Target  Models  #3  &  #4P(x)  =  4 e 

x2n  =~Y  (In  rn  +  In  rn+1)  =  average  cross 

section 


This  distribution  characterizes  a  target  with  one 
dominant  scatterer  plus  smaller  scatterers  or  one 
large  reflector  subject  to  small  changes  in  orien¬ 
tation  (Swerling  Cases  3  and  4) . 
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LIST  OF  RANDOM  NUMBERS 


The  following  list  of  random  numbers  is  suitable  for 
initializing  the  random  number  array  NRAND  used  by  the 
function  subprogram  RRAND.  The  numbers  were  derived  from 
tables  of  random  numbers  published  by  RAND  Corporation: 

A  Million  Random  Digits  with  100,000  Normal  Deviates, 
the  Free  Press  of  Glencoe,  New  York,  1955. 


*NL  1  0L  NR  A NO  =  1  2  06  3  1 5  8044 , 06  84  7664fc  5  9  ,  1  54  16  76  2 76 C  ,  19382  34  3 1  78  ♦ 
3  2122303427,  29052201E4C,  13*.  E 8647 C 55, Cl  734  7 3 1 40  E  , C 72 89355 507  , 

If  534467415 , 2438607  2834, 293 174535 72, C7 11464  364  2, 16232  718*23, 

CE  16C41  n’flC,  1641C917EI  3,2341652C791  ,2 d 82 36  38  45  2 , 10800745  445, 

29 M 76  16608, 2. 3 120  76  56  65,  32  320 4C266C, 194  71352757,0768375 551 7, 
C3f 65  7361  7  0  ,2 91  70  5  042  46 ,7 6 866 57  48  lb, 2 03 3 5bdC 6  12,  1486135  7  546, 

C 344 l  46  3H 2  2,2 5C 72 56 8248,2  1374  6 70C76 ,136)6667951,30463132192  , 
05804  7 76S 74  ,30172  034006,  If  18  4  26  1842,  149  74  2  1 64 67 , 1 C283C l  842C, 
302  56  5451  8b  ,1 33 1 0? 5 73 99 , 139 38 l b E? 07 , 0 12 b6C 746 5  7 , 196622 1 4 195 , 
103  32  795  36  1,015  7704  54  80,  1 6 74 2  £  6 165 5 , 1  lb E6 E 4 E  76 7  ,1 81 74 1 1 4680  , 
131741 14630 ,3089?48 7160, 3089248 7 160, 28360545700,33366415709, 
IC5567C  70 C  7,1  723592  1632  ,25322444850  ,30007056175  ,1348888  1553, 

10  140208896, 10224  1 4 £5  81, C  76554  2  3  38  7  ,2  26  26 40 2551 ,131010246  74  , 

05  7  79017119  ,30533512969  ,0721877  1539 ,00229 5 36t 7C , 25 1986Q44C 1 , 

Cl  70  2  666  304,151  34 181657, 27509(6 4464,203 2351 22 72  ,3  006822  7398  , 
170064588  73 ,1405100  739  3,  16 40 2  19 C290 , 2 f 3 C6 65 02 12 , 1 126C 71 7646 , 
168414 £27  74  ,1680162  6773 ,1  1  349273156 ,19404991345,069777128  30, 
2647  3264  7  2  1,0  288  34  3*1  3  7,  2  3  025  57CC91 , 1 10  1229 1622  ,13431365861  , 
17160  2929  17 ,31267410036, l 346 2 1 3 52 50 , 2646 3E 8 55 C 2 , 2 42 15 7 7 4 256 , 
25620  74  4156,1155782  06  55,3  0512  80  5  719  ,126  30  506319  ,17  722780814, 
03003554486, 0472255  '(022,  1650023  CO  91,16242624041, 16388044606, 

0  36659  53  72  3  ,2  62  12 6  9  84 3  8 , 1 3 5 70004  754 , 1  1 1 6  3 305  5 28 , 29 1 34 2 3 78 2 1 , 
16068  601392,1  2  l 64 54 2 C56 , 2 590 8 9 6 e 2 5 6  ,  C 356 4 S £66 66  ,245 1342 6529  , 
05383  8  738  59  ,25262  30  799  2,  16  4  It  25  17 7 7 , 3 2  1  4  5 2 705 39 , 2  11161 7  £ 576  , 
148247 11 U8C  ,193951  7304 3 , 2 0 74 306 11 71  ,2 13 19 355579 , 1907449  1967, 
180  81451743,  15244  3903  2  4, CEE 4f  12  235b  ,2  71 4c 30 5554  ,15825176938  , 
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SECTION  9 


RADS  IM  COMPUTER  PROGRAM 
FLOW  CHARTS 

PUBLISHED  IN  PART  2  OF  THIS  VOLUME. 


9-1  through  9-234 


SECTION  10 


INDEX 


(See  also  Section  6  in  Volume 

I,  Part  1 

and  Appendix  B  in  Volume  IV, 

Part  1) 

PROGRAM/ 

DESCRIPTION 

FLOW 

CROSS 

SUBROUTINE/ 

MODULE  &  LISTING 

CHART 

REF. 

FUNCTION 

NUMBERS  PAGE 

PAGE 

PAGE 

ABORT 

8-2 

9-81 

9-218 

ADDA 

207/239  See  CONV 

ADDEND 

235/236  See  RNDARY 

ADRNDC 

237 

See 

RNDARY 

ANTARY 

413 

Vol  III 

9-170 

9-230 

ANT TNT 

- 

8-4 

9-140 

9-226 

ANTPAT 

504 

See 

ANTARY 

ATOD 

216/217 

4-2 

9-98 

9-220 

AZGAIN 

- 

8-8 

9-159 

9-229 

BLOCK  DATA 

- 

8-10 

9-124 

9-224 

CDFNCL 

404 

See 

CDIGFL 

CDIGFL 

403 

4-6 

9-96 

9-220 

CFAR 

459,460 

4-10 

9-66 

9-216 

CGEN 

425 

Vol  III 

9-182 

9-231 

CGENCW 

455 

See 

CGENSF 

CGENSF 

453 

Vol  III 

9-149 

9-227 

CLINT 

302 

6-2 

9-198 

9-234 

CLUTTR 

503 

4-14 

9-112 

9-222 

CONS 

204 

5-2 

9-60 

9-216 

CONVMP 

205 

See 

CONV 

CUMDIS 

208/209 

6-11 

9-152 

9-228 

CUM  2 

- 

See 

CUMDIS 

DBLKX 

114/115 

8-12 

9-122 

9-223 

DCFAR 

440/441 

4-20 

9-93 

9-220 

DFT 

201 

4-24 

9-142 

9-226 

DFTF0 

234 

See 

DFT 

DFTRF 

233 

See 

DFT 

DIGFIL 

461 

4-30 

9-101 

9-221 

DIGFNC 

462 

See 

DIGFIL 

DIGFSF 

463 

4-34 

9-107 

9-222 

DIGTFL 

422/423 

4-41 

9-95 

9-220 

DIVA 

206 

See 

CONV 

. 

DTOA 

219/ 

220/228 

See 

ATOD 
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SUBROUTINE/ 

function 


ECM 

ECMFL 

ECMFSU 

ELGA IN 

ERGYCP 

ERGYRE 

ERRMSG 


MODULE 

NUMBERS 


512 


118 

116/117 
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&  LISTING 

PAGE 
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PAGE 
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REF. 

PAGE 

4-45 

9-104 

9-221 

DIGFIL 

DIGFSF 

8-8 

6-16 

9-160 

9-70 

9-229 

9-217 

See  ERGYCP 
See  ABORT 


FCENMP 
FGENXY 
F ILT 
FVDET 

HEX 

KLIM 

HWDET 

IBOOL 

IFLD 

TFT 

IFWDET 

IHLIM 

1  HWDET 
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LAMPCP 
LAMP  RE 
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MTIFLT 
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PACK 

PDF 


421 

420 

407/408 

416/417 

See 

See 

FGENXY 

4-51 

4-60 

HWDET 

9-72 

9-55 

9-217 

9-215 

4-64 

9-103 

9-221 

438/439 

414/415 

See 

HWDET 

4-65 

9-137 

9-225 

- 

8-14 

8-15 

9-119 

9-123 

9-223 

9-224 

See 

ZFFT 

447/448  See  HWDET 
442/443  See  HWDET 


445/446 

409/410 

511 

See 

HWDET 

4-70 

4-74 

8-16 

9-97 

9-117 

9-118 

9-220 

9-223 

9-223 

449/450 

See 

HWDET 

458 

4-80 

9-68 

9-217 

456/457 

See 

LAMPCP 

2-1 

9-2 

9-206 

2-2 

9-34 

9-210 

430/431 

4-83 

9-201 

9-234 

— 

See 

MTIFLT 

_ 

See 

MTIFLT 

432/433 

See 

MTIFLT 

436/437 

401/402 

See 

SWPINT 

4-91 

9-135 

9-225 

210/211 

See 

CUMDIS 

8-18 

9-88 

9-219 

212/213  See  CUMDIS 
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DESCRIPTION 
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FUNCTION 

PHDEC 

NUMBERS 

PAGE 

508/509 

4-94 

PHENC 

506/507 

4-100 

PLOTTR 

307-310 

6-18 

PLTFMT 

113 

6-23 

PTLIST 

303-306 

6-28 

PXFRM 

454 

Vol  III 

RDIGFL 

405 

4-105 

RECF 

451 

Vol  III 

RECFTF 

452 

See  RECF 

rndary 

214/215 

4-109 

RRAND 

101/102 

8-21 

RSHFTS 

231/232 

See  SHIFT 

RSHIFT 

229/230 

See  SHIFT 

RTOPDB 

103 

6-35 

RTOPM 

110 

See  RTOPDB 

RTOPM2 

111 

See  RTOPDB 

SCANNR 

313 

6-40 

SDBLKX 

- 

See  DBLKX 

SHIFT 

224/225 

4-113 

SHIFTS 

226/227 

See  SHIFT 

SPCAVG 

312 

Vol  III 

SQDET 

418/419 

See  HWDET 

SWPINT 

434/435 

4-118 

TARGET 

501 

4-122 

TGTNCL 

502 

See  TARGET 

TSARM 

- 

See  TSARY 

TSARY 

426 

Vol  III 

TSARY  1 

427 

See  TSARY 

TSARY  2 

- 

See  TSARY 

TSRPAT 

505 

Vol  III 

WEITCP 

222 

See  WEITRE 

WE  ITMP 

223 

See  WEITRE 

WEITRE 

221 

4-129 

WVGUID 

510 

4-135 

XYTODB 

104/108 

See  RTOPDB 

XYTOM 

105 

See  RTOPDB 

XYTOM2 

106 

See  RTOPDB 
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9-175 

9-82 

9-147 

9-100 

9-145 
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9-57 
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9-218 

9-229 
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9-231 

9-218 

9-227 

9-221 

9-227 

9-227 

9-223 


9-228 


9-229 

9-216 

9-229 

9-234 

9-222 


9-232 


9-231 


9-215 

9-223 
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PAGE 

ZFFT 

202 

5-7 

9-125 

9-224 

203  See  ZFFT 


SECTION  11 


RECURSIVE  COMPUTATION 
OF  SINE/COSINE  PAIRS 


Many  computer  programs  require  the  evaluation  of  sine 
and  cosine  functions  for  angles  that  are  uniformly  spaced 
over  the  interval  from  0°  to  360°,  for  example,  DFT  subroutine. 
In  many  cases  the  CPU  time  required  to  execute  these  com¬ 
puter  programs  can  be  significantly  reduced  by  the  use  of 
the  following  digital  oscillator  to  generate  sine/cosine 


pairs.  The  Z-plane  representation  of  the  digital  oscilla- 


Figure  11-1  Z-PLANE  REPRESENTATION  OF  THE  DIGITAL  OSCILLATOR 


The  variable  A 0  is  the  angular  increment  for  which  the 
sine/cosine  pairs  are  to  be  calculated.  The  digital  oscill¬ 
ator  block  diagram  is  shown  in  Figure  11-2. 

The  stimulus [cos  0O + j  sin  90]is  applied  only  when 
n=0  and  is  zero  otherwise.  Before  execution  of  the  digital 
oscillator  routine,  the  storage  registers  (A  &  B)  used  to 
generate  the  unit  delays  must  be  cleared. 


11-1 


Cos(nA0+9Q) 


Sin(nA9+0o) 


Figure  11-2  DIGITAL  OSCILLATOR  BLOCK  DIAGRAM 

n  =  number  of  executions  of  the  digital  oscillator 
routine . 

The  Fortran  statements  for  implementing  the  digital 
oscillator  are  the  following: 

A  =  Cos  (0O) 

B  =  Sin  (0O) 

DELC  =  Cos  (A 0) 

DELS  =  Sin  (AH) 


Loop 
n=  l ,  N 


A  and  B  contain  the 
sine /cosine  of  Oo-K 


n- 1  )Afl] ) 


TEMP  =  A 

A  =  A  *  DELC  -  B  *  DELS 
B  =  TEMP  •>  DELS  +  B  *  DELC 


where  A9  ~  nn.ule  increment. 

9,.  =  start  i  ng  angle  . 


MISSION 

of 

Rome  Air  Development  Center 


RAX  plans  and  conducts  research ,  exploratory  and  advanced 
development  programs  in  command,  control ,  and  communications 
(C3)  activities,  and  in  the  C3  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  communications,  electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation ,  solid  state  sciences ,  microeeve 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


